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In the baseline design for advanced LIGO interferometers, the most serious noise source is tiny, 
dynamically fluctuating bumps and valleys on the faces of the arm-cavity mirrors, caused by random 
flow of heat in the mirrors' sapphire substrates: so-called thermoelastic noise. We propose replacing 
the interferometers' baseline arm-cavity light beams, which have Gaussian-shaped intensity profiles 
that do not average very well over the dynamical bumps and valleys, by beams with mesa-shaped 
profiles that are flat in their central ~ 7 cm of radius, and that then fall toward zero as quickly 
as is allowed by diffraction in LIGO's 4 km arms; see Fig. |21 The mesa beams average the bumps 
and valleys much more effectively than the Gaussian beams. As a result, if the mirrors' substrate 
radii and thicknesses are held flxed at 15.7 cm and 13 cm, and the beam radii are adjusted so 
diffraction losses per bounce are about 10 ppm, replacing Gaussian beams by mesa beams reduces 
the thermoelastic noise power by about a factor 3. If other thermal noises are kept negligible, 
this reduction will permit advanced LIGO to beat the Standard Quantum Limit (circumvent the 
Heisenberg Uncertainty Principle for 40 kg mirrors) by about a factor 1.5 over a bandwidth about 
equal to frequency; optical (unifled quantum) noise will become the dominant noise source; and the 
event rate for inspiraling neutron star binaries will increase by about a factor 2.5. The desired mesa 
beams can be produced from input, Gaussian-profile laser light, by changing the shapes of the arm 
cavities' mirror faces from their baseline spherical shapes (with radii of curvature of order 60 km) to 
Mexican-Hat {MH; sombrero-like) shapes that have a shallow bump in the center but are otherwise 
much flatter in the central 10 cm than the spherical mirrors, and then flare upward strongly in 
the outer 6 cm, like a sombrero; Fig. |21 In this paper we describe mesa beams and MH mirrors 
mathematically and we report the results of extensive modeling calculations, which show that the 
mesa-beam interferometers are not substantially more sensitive than the baseline Gaussian-beam 
interferometers to errors in the mirror figures, positions, and orientations. This has motivated the 
LIGO Scientific Community (LSC) to adopt MH mirrors and mesa beams as an option for advanced 
LIGO, to be studied further. The details of our modeling calculations are presented in companion 
papers. 



I. INTRODUCTION AND SUMMARY 

The Laser Interferometer Gravitational- Wave Obser- 
vatory (LIGO) is designed to support successive gen- 
erations of interferometric gravitational-wave detectors. 
LIGO's first interferometers are now in operation 0, 
and the (negative) results of its first gravitational-wave 
searches have recently been published, j^]. When they 
reach their design sensitivity (presumably next year), 
LIGO's initial interferometers, together with their in- 
ternational partners, will reach out into the universe to 
distances where it is plausible, but not probable to de- 
tect gravitational waves Q. After a planned upgrade 
to advanced LIGO interferometers (planned to begin in 
2007), wave detection will be quite probable ^Sj. A base- 
line design for the advanced LIGO interferometers has 
recently been adopted along with several options, 
not currently in the baseline, that merit further study 
and might be incorporated at a future date. This pa- 
per describes one of these options, which has been much 
discussed within the LIGO Scientific Community (LSC) 
but has not previously been presented in the published 
literature: the reshaping of the arm-cavity light beams 
so as to reduce thermoelastic noise. 



A. The Context: Noise in Advanced LIGO 
Interferometers 

For advanced LIGO's baseline design the domi- 
nant noise sources in the most interesting frequency range 
(above about 20 Hz) are thermoelastic noise and optical 
noise (also called "unified quantum noise"). Other ther- 
mal noises (most especially coating thermal noise) might, 
in the end, be important; but in this paper we shall as- 
sume them negligible and shall focus on the thermoelastic 
noise and optical noise. 

In Fig. n we show the thermoelastic noise 0, @, the 
optical noise 0, 0, 0' their sum (labeled total noise), 
and the standard quantum limit (SQL) for the advanced 
LIGO baseline design with sapphire mirrors This 
figure suggests (as is well known J\) that, if the ther- 
moelastic noise can be reduced significantly (and if other 
thermal noises can be kept negligible), then the advanced 
LIGO interferometers will be able to beat the SQL, and 
the interferometers' ranges (detectable distances) for as- 
trophysical sources will be increased significantly. 

In this paper we propose a method ("fiattening the 
interferometers' light beams") for reducing the thermoe- 
lastic noise, we evaluate the resulting increased range for 
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FIG. 1: Noise curves for advanced LIGO. 



neutron-star / neutron-star (NS/NS) binaries, and we ex- 
plore practical issues related to our proposal. We have 
previously discussed our proposal, the increased NS/NS 
range, and the practical issues in presentations at meet- 
ings of the LIGO Scientific Collaboration [Hill [Hill 
and in an internal LIGO document 

It is worth emphasizing that our proposed beam flat- 
tening need not entail any major modifications of the 
advanced-LIGO interferometer design; it only involves 
reshaping the reflecting surfaces of the test mass mirrors, 
adjusting the input beam, and modifying the mirror- 
alignment control system to account for the changes of 
the hght-mode shapes. 



B. The Physical Nature of Thermoelastic Noise; 
Motivation for Reshaping Beams 

Our proposal is motivated by the physical nature of 
thermoelastic noise. This noise is created by the stochas- 
tic flow of heat (random motions of thermal phonons) 
within each test mass (mirror), which produces stochas- 
tically fluctuating hot spots and cold spots inside the test 
mass. The test-mass material (sapphire for the baseline 
design of advanced LIGO) expands in the hot spots and 
contracts in the cold spots, creating fluctuating bumps 
and valleys on the test-mass (mirror) faces. These face 
bumps influence the light beam's measurement of the test 
masses' positions: the interferometer's output phase shift 
is proportional to the difference of the test masses' aver- 
age positions — with the average being the position of 
a mirrored test-mass face, weighted by the light's energy 
flux (its intensity distribution). 

If the intensity distribution is "flat" (nearly constant) 
in most regions of high intensity, then the adjacent valleys 
and bumps (having been created by heat flow from one to 
the other) will average out, giving low net thermoelastic 



noise. If, instead, the energy flux is sharply changing in 
most regions of high flux, then the adjacent valleys and 
bumps will not average well and the thermoelastic noise 
will be high. Also, the larger is the light beam, the better 
will be the averaging and thus the lower will be the noise. 

These considerations suggest that large-radius, flat- 
topped beams with steep edges (e.g. the thick curve in 
Fig. 121 below) will lead to much smaller thermoelastic 
noise than small-radius, centrally peaked beams with 
gradually sloping sides (e.g., the thin, Gaussian curve 
in Fig. H below). 

For Gaussian beams, the influence of beam radius Tq 
has been quantified by Braginsky, Gorodetsky and Vy- 
atchanin j^] (who first pointed out the importance of 
thermoelastic noise for sapphire test masses): the ther- 
moelastic noise power scales as S"™ (x l/rj^ (aside from 
small corrections due to the test masses' finite sizes 0). 
This has motivated the baseline design for advanced 
LIGO interferometers with sapphire test masses: the 
beam radius ro is chosen as large as possible, given the 
demand for small diffraction losses, £o ^ 10 ppm per 
bounce in the interferometer's arm cavities.^ 

The baseline design uses light beams with a Gaussian 
distribution of energy flux, since such beams are eigen- 
functions of cavities with spherical mirrors, and spherical 
mirrors are a standard, well-developed technology. How- 
ever, the Gaussian energy flux is far from flat: most of 
the energy is in regions where the flux is rapidly varying 
with radius (thin curve in Fig. |21 below), and correspond- 
ingly the thermoelastic noise is substantially larger than 
it would be with "flat-topped" beams (thick curve in Fig. 
12}. This has motivated a (previously unpublished) pro- 
posal by O'Shaughnessy and Thorne 10] to replace the 
Gaussian beams with flat-topped beams, while keeping 
the beam radius as large as is compatible with diffrac- 
tion losses Cq < 10 ppm. 



C. Summary of Analysis and Results 

In Sec. m] we construct an example of a flat-topped 
light beam — a flat-topped TEMOO mode of light that 
will resonate in an interferometer's arm cavity, if the test- 
mass mirror faces are shaped appropriately. Because the 
intensity distribution of our flat-topped beam resembles a 
mesa in the southwest American desert, we call it a mesa 
beam^ (a name suggested to us by Phil Willems). To 
produce this mesa beam as an eigenmode of a symmetric 



^ The 10 ppm is dictated by the following considerations: For the 
baseline design there is 125 W of input power to the interferom- 
eter and 830 kW of circulating power in each arm cavity. Ten 
ppm of diffraction loss per bounce results in a diffraction power 
loss in the arm cavities of 4 X lOppm X 830kW = 33 W, which is 
25 per cent of the 125 W of input light, a reasonable value. 

^ It is also called a flat-topped beam and a mexican - hat or MH 
beam in the internal LIGO literature [Tol [Tl [TllTl [ll 
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arm cavity one must give the mirror faces a shape, with 
a central bump and an upturned brim, that resembles a 
Mexican hat (or sombrero) (Fig.OJ, so we call the mirrors 
Mexican-hat (MH) mirrors. 

We have not optimized our mesa beams' intensity dis- 
tribution so as to bring the thermoelastic noise to the 
lowest value possible, but in Sec. ^] we argue that our 
chosen mesa beams are likely to be close to optimal. 

Figure|21below shows the intensity distribution for our 
proposed mesa beam (thick curve) and compares it with 
the intensity distribution of the baseline Gaussian beam 
(thin curve), which has the same diffraction losses. Fig- 
ure 13 below compares the mirror shapes that support 
these mesa and Gaussian beams as eigenmodes of a 4 km 
LIGO arm. In their inner 10 centimeters of radius, the 
MH mirrors that support mesa beams are much flatter 
than the spherical mirrors that support Gaussian beams, 
but in their outer 6 centimeters (the upturned brim re- 
gion), the MH mirrors arc far more curved. 

As we shall see, this greater curvature at large radii 
compensates considerably for the flatter shape at small 
radii, enabling a mesa-beam interferometer to exhibit 
only modestly worse parasitic-mode behavior than a 
Gaussian-beam interferometer, and only modestly worse 
sensitivity to mirror tilts, displacements, and figure er- 
rors. 

Three of us (O'Shaughnessy, Strigin and Vyatchanin; 
OSV 15]) have computed the substantial reductions in 
thermoelastic noise that can be achieved in advanced 
LIGO by replacing the baseline spherical mirrors and 
their Gaussian beams with MH mirrors and their mesa 
beams. The method of computation and the results are 
described in Sec. IIIII Our principal conclusion is this: 

1. By switching from the baseline (BL) spherical mir- 
rors to MH mirrors with the same cylindrical test- 
mass diameters and thicknesses and the same 10- 
ppm-per-bounce diffraction losses as the BL, one 
can reduce the power spectral density of thermoe- 
lastic noise by a factor 0.34 and increase the volume 
of the universe reachable for compact-binary inspi- 
rals, and thus also LIGO's binary inspiral event 
rate, by a factor 2.6. Larger improvements could 
be achieved by using conical test masses with en- 
larged inner faces (frustums; Tables HI and IV and 
Appendix I of Xt'D- 

One might worry that the greater flatness of the MH 
mirrors, in the inner 10 cm where most of the light re- 
sides, will make mesa-beam interferometers much more 
sensitive to errors in the orientations, positions, and fig- 
ures of the mirrors. We have explored this issue in great 
depth, with the conclusion that mesa-beam interferome- 
ters are not substantially more sensitive to mirror errors 
than Gaussian-beam interferometers. Details of our ex- 
plorations are given in companion papers by D'Ambrosio 
|l6j | and by O'Shaughnessy, Strigin and Vyatchanin [Tsj . 
and our methods and conclusions are presented and dis- 
cussed in Sec. IIVI of this paper. Our quantitative con- 



clusions rely on a comparison of two fiducial advanced- 
LIGO interferometer configurations, both of which have 
the same L — 4 km arm lengths, and the same R ~ 16 cm 
coated mirror radii, and beam radii that give the same 
18ppm diffraction losses off the mirrors. One configura- 
tion uses fiducial spherical mirrors with radii of curva- 
ture TZc = 83 km and so has g = 1 — L/TZc = 0.952, 
and supports fiducial Gaussian beams with beam radii 
at the mirrors = 4.70 cm. (This configuration differs 
from the current baseline design for advanced LIGO; it 
is nearer an altered baseline that we recommend in Sec. 
nil A 21 ) Our second configuration uses fiducial MH mir- 
rors and fiducial mesa beams with beam radii (defined in 
Sec. Ill A|) D — ib — 10.41 cm. Our comparative conclu- 
sions for these two configurations are the following (see 
Sec. IIVI for greater detail). [Note: if we had compared 
the baseline Gaussian-beam configuration with a mesa- 
beam configuration with the same diffraction losses, we 
would have reached approximately the same conclusions; 
see, e.g., the paragraph following Eq. 

2. Among those parasitic optical modes of a perfect 
arm cavity that are not strongly damped by diffrac- 
tion losses, the parasite closest in frequency to the 
desired TEMOO mode is separated from it by 0.099 
of the free spectral range in the fiducial Gaussian- 
beam case, and by 0.0404 of the free spectral range 
in the fiducial mesa-beam case; see Sec. lIVd This 
factor ^ 2 smaller mode spacing leads to a modestly 
greater sensitivity of the mesa-beam interferometer 
to mirror tilt errors. 

3. The interferometer's arm cavities are about four 
times more sensitive to mirror tilt when fiducial 
MH mirrors are used than for fiducial spherical 
mirrors. When all four cavity mirrors are tilted 
through angles 9 about uncorrelated axes, the frac- 
tions of the carrier power driven into (dipolar) par- 
asitic modes inside the arm cavities, and driven 
out the dark port, are about 0.001(61/0.01 /xrad)^ 
and 0.002(61/0.01 ^irad)^ respectively for our fidu- 
cial MH mirrors; and 0.002(6/0.035 /irad)^ for our 
fiducial spherical mirrors; Sec. IIVDI This factor 
four greater sensitivity is not a serious issue, since 
it turns out that the strictest constraints on mirror 
tilt come from the signal recycling cavity (and, if a 
heterodyne output were to be used, from the power 
recycling cavity), and not from the arm cavities; see 
items [71 and |H1 below. 

4. The interferometer's arm cavities are about equally 
sensitive to transverse displacements of their end 
test-mass mirrors (ETM's), whether the mirror 
shapes are fiducial MH, or are fiducial spherical. 
Specifically: for uncorrelated displacements of the 
two ETMs through distances s, the fractions of the 
carrier power driven into (dipolar) parasitic modes 
inside the arm cavities, and driven out the inter- 
ferometer's dark port, are about 100(s/1.0mm)^ 
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ppm and 190(s/1.0nim)^ ppm, respectively, for 
fiducial MH mirrors; and 100(s/1.3mni)^ Ppm and 
190(s/1.3mm)^ ppm for fiducial spherical mirrors. 
For details, see Sec. IIVEI 

5. For MH mirror figure errors with peak-to- valley 
height variations Az in the innermost 10 cm by 
radius: after the control system has optimized 
the mirror tilts, the fractions of the carrier power 
driven into parasitic modes inside the arm cav- 
ities, and driven out the dark port, are about 
0.0008(Az/6nm)2 and 0.0015(Az/6 nm)2, respec- 
tively; Sec. IIV Fl For spherical mirrors, a figure er- 
ror Az about twice as large drives the same power 
into parasitic modes — and the relevant figure error 
is confined to a smaller central region of the mirror 
(about 7.3 cm radius for baseline spherical mirrors 
and 8 cm for fiducial spherical mirrors, compared 
to 10 cm for fiducial MH mirrors). The measured 
mirror figure errors in the initial LIGO interferome- 
ters are of order Az = 3 to 6 nm peak to valley (i.e. 
1 to 2 nm rms), which suggests that the MH arm 
cavities' required figure errors may be achievable. 

6. The most serious constraints on mirror tilt and on 
mirror figure accuracy come not from the arm cav- 
ities but rather from the signal recycling (SR) cav- 
ity. The SR cavity, power recycling (PR) cavity, 
and RF sideband cavity operate approximately in 
the geometric optics regime and thus are nearly in- 
sensitive to whether one uses MH or spherical mir- 
rors; Sec. IIV HI As a result, by switching from 
spherical to MH mirrors, one pays only a small 
penalty, in terms of mirror tilt constraints and 
figure-error constraints. 

7. More specifically, the most severe constraints on 
tilt and figure error arise from the driving of signal 
power into parasitic modes when the signal light 
passes through the SR cavity. To keep the re- 
sulting loss of signal strength below one per cent 
in the standard wideband advanced LIGO inter- 
ferometers, it is necessary to constrain the magni- 
tude 6 of the vectorial tilts of the input test-mass 
mirrors (ITM's) and signal recycling mirror (SRM) 
to 61^3 < 0.024 /xrad (for the fiducial spherical 
mirrors) and 6*^3 < 0.016 /irad (for the fiducial 
MH mirrors). When an advanced interferometer 
is narrowbanded at / ~ 500 Hz or ~ 1000 Hz, 
the constraint must be tighter: 6*^3 < 0.011 fiiad, 
and < 0.007 /xrad. These are approximately 
the same as the experimentally achieved ^3 con- 
straints on LIGO-I tilt arising from the RF cavity 
in the absence of an output mode cleaner, and the 
same as the baseline constraints on tilt in advanced 
LIGO arising from misalignment of the input laser 
beam[l8l|. If there were no output mode cleaner 
in advanced LIGO and heterodyne readout were 
used in place of the baseline homodyne readout. 



then the constraint on tilts in the RF cavity (due 
to mode mixing for the RF sideband light used in 
the readout) would be about the same as that for 
narrowband interferometers in the SR cavity. For 
the baseline homodyne readout, no such RF con- 
straint arises. The reduction in signal strength due 
to all the above (item[7|l tilts scales as 9^; and we 
estimate that the above (item [7|l constraints are 
inaccurate by a factor < 2 due to ignoring corre- 
lations in the overlaps of certain parasitic modes, 
and for the narrowbanded interferometers, due to 
inaccuracy of the geometric optics approximation 
in the SR cavity. For details of all these issues, see 
Sec. HVTI 

8. We characterize the analogous constraints on mir- 
ror figure error by the peak-to-valley fluctuations 
in the mirror height in the central regions of the 
mirrors (regions enclosing 95 per cent of the light 
power; radius ~ 10 cm for fiducial MH mirrors and 
~ 8 cm for fiducial spherical mirrors), with the 
fluctuations averaged over ~ 3 cm (an averaging 
produced by breakdown of geometric optics in the 
SR cavity). Our estimated constraints for one per 
cent loss of signal strength are Azwb ^ 2.0 nm 
for wideband advanced LIGO interferometers and 
Aznb ^1-0 for narrowband, independently of 
whether the mirrors are MH or spherical — though 
the region over which the constraints must be ap- 
plied is different, 10 cm radius for MH and 8 cm for 
spherical. The loss of signal strength scales as Az^, 
and our estimated constraints might be inaccurate 
by as much as a factor ~ 3 due to exploring only 
one representative shape for the figure errors, due 
to overlaps of certain parasitic modes, and for the 
narrowbanded interferometer due to inaccuracy of 
the geometric optics approximation in the SR cav- 
ity. These are approximately the same constraints 
as arise (in our calculations) from the RF cavity in 
LIGO-I, in the absence of an output mode cleaner. 
If there were no output mode cleaner in advanced 
LIGO and heterodyne readout were used, then the 
constraint on figure errors in the RF cavity (due to 
mode mixing for the RF sideband light used in the 
readout) would be about the same as that for wide- 
band interferometers in the SR cavity. For details 
of these conclusions, see Sec. IIV Jl 

Among all the constraints on mirror errors that arise 
from our modeling, the most serious are the last ones 
(items and |SJ : SR-cavity- induced constraints on mir- 
ror figure errors to avoid a one per cent loss of signal 
strength. These constraints are nearly independent of 
whether the mirrors are spherical or MH. These con- 
straints would be relaxed if the SR cavity were made less 
degenerate. This could be achieved by shaping the fronts 
of the ITMs as lenses that bring the light (Gaussian or 
mesa) to a focus somewhere near the SR mirror — and 
also near the PR mirror. 
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Because MH mirrors and their mesa beams produce 
such a great (factor 3) reduction of thermoelastic noise 
power, and they increase the sensitivity to mirror errors 
by only modest amounts, they have been adopted as op- 
tions for advanced LIGO, and they may be of value for 
LIGO's future international-partner interferometers. In 
Sec. we describe some of the future research that is 
needed in order to firm up our understanding of the pros 
and cons of MH mirrors and mesa beams. 



D. Notation 

We here summarize some of the notation used in the 
remainder of this paper. The numerical values are for 
advanced LIGO interferometers, including sapphire test- 
mass substrates, with the sapphire idealized as isotropic 
(its properties averaged over directions).^ 

b: Diffraction lengthscale b — ^ \L/2'k — 2.603 cm for 
light with wavelength A = 1.064 /im in the L = 4km 
LIGO beam tubes; equal to a symmetric Gaussian 
beam's minimum possible radius at the end mirrors 
(the radius at which the power flux has dropped to 
1/e of its central value). 

Cv '• Specific heat of test-mass substrate per unit mass at 
constant volume [7.9 x lO^erg g^^ K^^] 

D: Mesa beam radius [units cm]; Eq. (|2.3|l 

E: Young's modulus of test-mass substrate [4 x 
10^^ dynccm^^] 

ETM: End test mass of an arm cavity 

H: Thickness of test mass [units cm] 

/: Gravitational-wave frequency at which noise is evalu- 
ated [units Hz] 

g: Arm cavity's g-factor, g = \ — L/TZc where L is arm 
length and TZc is the (identical) radius of curvature 
of its mirrors 

J-: Finesse of an optical cavity 

/: Noise integral for a test mass, Eq. H3.2() 

ITM: Input test mass of an arm cavity 

k: Wave number, equal to 2ti/\ 



^ The influence of mechanical anisotropies and of heat-conduction 
anisotropics on the thermoelastic noise have not yet been inves- 
tigated, but could be large enough to be important, since the rel- 
evant parameters vary with direction by amounts of order ±5% 
for the Young's modulus E, ±25% for 1 — 2a, ±5% for the heat 
conductivity k, and ±5% for the thermal expansion coefficient 
ai. 



ks- Boltzmann's constant [1.38 x 10"-^^ ergK""^] 

L: Interferometer arm length [4 x 10''^cm] 

C: Diffraction loss in a single reflection off a mirror 

M : Mass of test mass [4 x lO^'g] 

PRM: Power recycling mirror 

r Radius in transverse plane 

Tq: Radius, on test-mass face, at which the intensity of 
a Gaussian light beam has dropped by a factor 1/e 
from its central value [ro = 4.23 cm = 1.636 for 
baseline design]. Note: many LIGO papers, fol- 
lowing SiegmannfT^. define the beam radius to be 
V2 larger, so it is the radius at which the beam's 
amplitude has dropped by 1/e. 

P„: Fraction of interferometer's light power in mode n 

Rpi The physical radius of a test mass [15.7cm] 

R: The radius of the mirror coated onto a test mass 
[equal to Rp or Rp — 8mm]; also, the power re- 
flectivity of a mirror 

TZc Radius of curvature of a spherical mirror [units cm] 

s: Transverse displacement of an arm cavity's ETM 

Sh{f)- Spectral density of noise (thermoelastic or other) 
for detecting a gravitational wave h with optimal 
direction and polarization [units Hz"-'^] 

SRM: Signal recycling mirror 

T: Temperature of test-mass substrate [300K] 

u: Electric field of some light mode or superposition of 
modes (renormalized to unit norm, J ]u]^c?Area= 
1); usually evaluated at the transverse plane tan- 
gent to an ITM mirror face, with the light propa- 
gating away from the ITM. Subscripts identify the 
mode. 

U: Unnormalized electric field of some light mode. 

v: Same as u: Unit-normed electric field of some light 
mode or superposition of modes. 

Az: The peak-to-valley mirror deformation (mirror fig- 
ure error) in a mirror's central region 

ar. Substrate's coefficient of linear thermal expansion 
[5.0 X 10"*^ K-i] 

ai^2'- Amplitude of excitation of an arm cavity's parasitic 
mode ui,2 by a tilt of the cavity's ETM; Eq. H4.3|) 

/3i : Amplitude of excitation of the parasitic mode vi by 
mirror figure errors; Eq. H4.16|l 



6 



7o Overlap of arm cavity's fundamental mode uq with 
Gaussian mode Ud that drives it; Eq. 14. 2|) 

6i Fraction of the light power of some perturbed field u'^ 
that is in parasitic modes; Eq. (|4.29() 

k: Thermal conductivity of test-mass substrate [3.3 x 10^ 

erg s^^ cm^^ K~^] 

A: Wavelength of laser light [1.064 /im]; also, in Sec. 
II V 1 31 a function appearing in the analysis of the 
signal recycling cavity. 

p: Density of test-mass substrate [4gcm~^]; also, ampli- 
tude reflectivity of signal recycling mirror 

a: Poisson ratio of test-mass substrate [0.23] 

6: Angle of mirror tilt 

6'diff: Diffraction angle of a Gaussian beam (the Gaus- 
sian beam's opening angle if it were to continue on- 
ward through the mirror to large distance); 6'diff ~ 
X/nWo = 2[(1 - g)/{l + g)]^/'^b/L, where Wo is the 
mode's waist radius (radius at which the amplitude 
has dropped by 1/e from its central value). 

Q: expansion (fractional volume change) of substrate 

Lu = 271 f: Angular frequency, corresponding to the fre- 
quency / at which the noise St is evaluated 

C,iX- Amplitude of excitation of an arm cavity's parasitic 
mode by transverse displacement of the ETM; Eq. 

(ion 



II. MEXICAN-HAT MIRRORS AND THE MESA 
MODES THEY SUPPORT 

In this paper we study a specific variant of a mesa light 
beam and the MH mirrors that support it. We believe 
this variant to be near optimal for reduction of thermoe- 
lastic noise, but we have not carried out the (rather com- 
plex) analysis required to prove optimality. 



A. Mesa Fields 

The flat-topped (mesa-shaped) eigenmode of an inter- 
ferometer arm cavity, which we seek to construct, must 
have an intensity distribution that is nearly flat across 
most of the light beam, and that then falls as rapidly as 
possible (constrained by diffraction effects) at the beam's 
edges. Moreover, if (as in baseline advanced LIGO) the 
cavity's input test mass (ITM) and end test mass (ETM) 
have the same physical dimensions, then to minimize the 
thermoelastic noise at fixed net diffraction loss, the beam 
should be symmetric about the arm cavity's mid plane, 
so its beam radii D are the same on the two mirrors. Oth- 
erwise [since 5™ oc 1 / approximately, and diffraction 



losses increase exponentially rapidly with increasing _D; 
Eq. the mirror with the reduced beam radius and 
smaller diffraction loss will have its thermoelastic noise 
power increased, while that with the enlarged beam ra- 
dius and larger diffraction loss will have its noise power 
decreased more modestly, leading to a net noise increase. 

The fastest possible falloff, for light in an optical cavity 
of length L, is that on the edge of the minimal Gaussian 
beam — the Gaussian beam whose radius increases by a 
factor y/2 in going from the beam waist (at the cavity's 
center plane) to the cavity's end mirrors. This minimal 
Gaussian, at the mirror planes, has the following (unnor- 
malized) form 



t^min Gauss =exp 



-r2(l + z) 



2&2 



where 



^/\L/2n = 2.603 cm , 



(2.1) 



(2.2) 



with L = A km the cavity length, k — 27r/A the wave 
number, and A = 1.064/im the wavelength of the light. 

A near-optimal flat-topped eigenmode, with near- 
minimal thermoelastic noise, will have a constant inten- 
sity in the central region, and will fall off at its edge at 
approximately the same rate as the edge of this minimal 
Gaussian. To produce such an eigenmode, we superpose 
minimal-Gaussian fields, with a field density that is con- 
stant out to a radius r = D and then stops abruptly. 
More specifically, our chosen unnormalized eigenmode 
has the following form: 



U{D,r) = 



exp 



2&2 



dxodyo 



(2.3) 



where r = ^ -V y"^ and the integration is over a circle 
Cu of radius D: < I?. 

By carrying out the yo integral in Cartesian coordi- 
nates, with 2/ = and a: = r, we obtain the following 
expression for U{D,r), which wc have used in much of 
our numerical work: 



U{D,r) 



-2n 
TTi 



dxo exp 



-{xo-rf{l + i) 



262 



xerfi 



-(1 + *) 



(2.4) 



Here erfi(z) — eri{iz)/i is the imaginary error function. 

By converting to circular polar coordinates and per- 
forming the angular integral, we obtain the following sim- 
pler expression for our unnormalized eigenmode 



U{D,r) 



2tt I exp 
'o 



-(r2+r2)(l + *) 



X/n 



rro{l + i) 
62 



262 

TodVo ■ 



(2.5) 
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Here Iq is the modified Bessel function of order zero. 
Modes with other weightings of the minimal-radius Gaus- 
sians can be obtained by inserting a weighting function 
/(fo) into the integrands of Eqs. (|2.3|l and (|2.5I) . 

In the Appendix we give some approximate formulae 
for U{D,r) valid at large radii. These are useful for 
quick, clipping-approximation computations of diffrac- 
tion losses. 

The squared norm of U{D, r) (the area integral of its 
squared modulus) is given by the following approximate 
formula, which is accurate to within a fraction of a per 
cent for 3.0 < D /b < 6.0 (the regime of interest to us): 

/•oc 

N^{D)= / |f/(i:),r)|227rrdr = 4.66-50. 58i:)-h62.10i:'^ 
Jo 

(2.6) 

We denote by u the normalized field on the mirror faces, 
and to distinguish it from a Gaussian field, we sometimes 
will use a subscript "mesa" : 



,(D,r) =u(Ar) 



NiD) 



(2.7) 



B. Gaussian fields 

The advanced LIGO baseline design uses arm cavi- 
ties with spherical mirrors, which have Gaussian modes 
whose field at the mirror plane is (cf. (19|) 



UGir,ro) = 



1 



: CXp 




(2.8) 

Here Tq is the beam radius (at which the energy flux falls 
to 1/e of its central value). From the phase of this field 
one can read off the radius of curvature of the mirrors 



7^c - i 



(2.9) 



and thence the arm cavity's g-factor 
g - 1 - L/Uc . 

C. Diffraction Losses 

In the baseline design of an advanced LIGO interfer- 
ometer 0, the test masses are cylinders whose faces are 
coated with dielectric mirrors out to a radius R that is 
0.8cm less than the cylinders' physical radii 



R = Rp- 0.8cm 



(2.10) 



We shall explore MH mirrors that are coated in this same 
manner, R = Rp — 0.8 cm and also MH mirrors that are 
coated all the way out to the test-mass edges, R — Rp. 
The diffraction losses in each reflection of a cavity mode 



off a mirror are given, approximately, by the clipping 
approximation 



POO 

Cciip = / \u{r)f2TTrdr 

J r 



(2.11) 



Here u(r) is the normalized field [Minosa(^, '') for a mesa 
mode and UG{i'o,'r) for a Gaussian mode with infinite 
mirrors, R — oo\. 

In actuality, the mirrors' edges at r = i? modify the 
field thereby causing the true diffraction losses to differ 
from this clipping formula. The true diffraction losses 
have been computed by OSV via a numerical solu- 
tion of the eigenequation for the cavity modes, and in- 
dependently by D'Ambrosio [T^ using an FFT code to 
propagate light in the cavity. The results are 

£o — 0.85£ciip for mesa modes , 

£o — 2.5£ciip for Gaussian modes , (2.12) 

in the parameter regime of interest — though the numer- 
ical coefficients 0.85 and 2.5 can oscillate substantially as 
the beam radii and mirror radii are changed. When we 
need high-accuracy diffraction losses (e.g. in portions of 
Sec. II V|) . we compute them with care using the cavity 
eigenequation 0) or FFT code 



D. Mirrors and Normalized Flux for Mesa and 
Gaussian Modes 



The baseline design for advanced LIGO interferometers 
has mirror radii R = 14.9 cm and Gaussian beam-spot 
radii rg = 1.636 = 4.23 cm, corresponding to a diffraction 
loss of £o = 10 ppm and a mirror radius of curvature 
TZc = 53.7 km. The normalized energy flux \ug{to,t)\'^ 
for this baseline Gaussian field is shown in Fig. [3 and 
the shape of the mirror (segment of a sphere with radius 
53.7 km) is shown in Fig. |31 

A cavity made from MH mirrors with the baseline 
radius R = 14.9 cm and the baseline diffraction losses 
Cq = 10 ppm has a mesa beam radius D — 3.435 = 8.92 
cm [computed from Eqs. H2.12|l and H2.11|l ]. The normal- 
ized energy flux |Minesa(£', t)\'^ for this mesa field is shown 
in Fig. 121 Notice how flat the top of this intensity proflle 
is, and how much like a mesa the profile is shaped, and 
notice the contrast with the Gaussian profile. 

The surfaces of the MH mirrors coincide with the mesa 
field's surfaces of constant phase; i.e., their height 5z as 
a function of radius r is given by 



k5z = Arg[umosa(£', r)] , 



(2.13) 



where k — 2n/X is the light's wave number. This MH 
mirror shape is shown in Fig.|3| Notice the shallow bump 
in the middle and the flaring outer edges. This bump and 
flare resemble a Mexican hat (sombrero) and give the MH 
mirror its name. 
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radius, cm 



FIG. 2: The power distributions of: (i) the Gaussian mode for 
the baseline mirrors with coated mirror radius R — 15.7 cm 
and beam radius Vo — 1.73b — 4.50 cm (thin curve), which 
has diffraction loss per bounce £o = 10 ppm; and (ii) the 
mesa mode with D = 3.73fe — 9.71 cm (thick curve) which, 
for this same coated mirror radius R — 15.7 cm, has the same 
diffraction loss per bounce £o = 10 ppm. 




2.5 5 7.5 10 12.5 15 
radius, cm 

FIG. 3: The shape of the spherical mirrors (thin curve) and 
MH mirrors (thick curve) that support the arm-cavity Gaus- 
sian mode (thin curve) and mesa mode (thick curve) of Fig.|21 
The height is measured in units of the wavelength of the light, 
A = 1.064/im. 



III. THERMOELASTIC NOISE AND 
NEUTRON-STAR BINARY RANGE FOR 
MESA-BEAM INTERFEROMETERS 

A. Thermoelastic Noise 

1. Quantifying the thermoelastic noise: the noise integral 
and the NS/NS range 

Building on the seminal work of Braginsky, Gorodetsky 
and Vyatchanin Liu and Thorne 6] have used Levin's 
[20(| direct method to derive the following formula for an 
interferometers' thermoelastic noise in terms of a noise 



integral I a- 

where / is the average, over the four test masses, of the 
thermoelastic noise integral, 

/ = i ^ , Ia = ^ f (Ve)2dvolume ; (3.2) 

A=l ° 

cf. Eqs. (3), (4) and (13) of In Eq. igU, the notation 
is as spelled out in Sec. II Dl and we use numerical values 
(shown in Sec. II Dl) that assume the test-mass substrate 
is sapphire, idealized as an isotropic material.l^ In Eq. 
(|3.2|) . 8 is the expansion (fractional volume change) in- 
side the test-mass substrate, produced by a static force 
with magnitude Fg and with profile identical to that of 
the light beam's intensity distribution over the test-mass 
face (e.g., Fig. [2Jl, and the integral is over the volume 
Va of test-mass A. Note that the dimensions of Ia and 
thence of / are length/force^ — s'*g~^cm^^. 

Equation 1)3.1(1 shows that the frequency dependence 
of the thermoelastic noise is independent of the mirror 
shape and test-mass shape; it always has the same slope 
as the SQL (except in testbed systems with tiny mirrors 
and light beams j^], which are irrelevant in this paper). 
As a result, the thermoelastic noise produced by an ad- 
vanced LIGO interferometer whose mirrors have some 
chosen shapes, divided by the thermoelastic noise of the 
baseline advanced LIGO interferometer, is equal to the 
ratio of the two interferometers' noise integrals 

[Shif)/Sf'^if)]TE = I/Ibl (3.3) 

This motivates our use of I / Ibl as one measure of a 
candidate interferometer's thermoelastic noise. 

We shall also use a second measure: The distance 
to which the LIGO network can detect NS/NS bina- 
ries, with network amplitude signal-to-noise ratio 8; we 
call this the network's range for NS/NS inspirals. This 
three-interferometer network range is larger by a factor 
-s/S = 1.732 than the single-4km-interferometer NS/NS 
range that is often used by the LIGO community and 
that is encoded into the "BENCH" LIGO software 
In computing the network range for advanced LIGO with 
Gaussian or mesa beams, we assume that; (i) the three 
advanced LIGO interferometers (all with L — A km) are 
all being operated with signal-recycling-mirror parame- 
ters that maximize the range (the operation mode ten- 
tatively planned for the first year of advanced LIGO ob- 
servations), (ii) they all incorporate identical sapphire 
mirrors with the same shapes and beam sizes, (iii) all 
thermal noises are negligible except thermoelastic noise, 
and (iv) the remaining interferometerparameters have 
their baseline advanced LIGO values 0| (e.g., the circu- 
lating power in each arm is 830 kW). For the baseline 
advanced LIGO design with Gaussian beams, the single- 
4-km-interferometer NS/NS range, computed under these 
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0.2 0.4 0.6 0.8 1 



(Sjj/Sj^ )te 

FIG. 4: The LIGO network NS/NS range as a function of 
the thermoelastic noise power, in units of the baseline ther- 
moelastic noise, {Sh/Sf^)TE- For each thermoelastic noise 
level, the advanced LIGO interferometer's optical parameters 
(homodyne readout phase and signal-recycling mirror) are op- 
timized to produce the greatest possible NS/NS range. The 
optimization has been performed for us by A. Buonanno and 
Y. Chen (private communication), assuming that the only 
significant noise sources are thermoelastic noise and optical 
(unified quantum) noise. 



2. Baseline test masses with spherical mirrors and 
Gaussian beams 

The baseline design of an advanced LIGO interferom- 
eter entails four identical test masses: sapphire cylinders 
with physical radii Rp = 15.7 cm, coated-mirror radii 
R ~ Rp — 0.8 cm, thickness H — 13 cm, density p = A 
g cm^"^ and mass M — AO kg; and the baseline light 
beam at the test-mass face is Gaussian with beam radius 
To = 4.23 cm = 1.636 so the diffraction loss per bounce 
is Co — 10 ppm. For this baseline beam and test mass, 
OSV ^3 find for the value of the noise integral 

/bl = 2.57 X lO^^^s^g^^cm^i . (3.5) 

We advocate extending the mirror coating out to the 
test- mass edge so R = Rp = 15.7 cm, and increasing the 
Gaussian beam radius correspondingly, to ro = 4.49cm — 
1.7256, so the diffraction losses are still 10 ppm. With 
this beam expansion, OSV find (Table I of 15]) that the 
thermoelastic noise is reduced to Sh/Sj^^ ~ I/Ibl = 
0.856 15], from which we deduce (via Fig. that the 
range for NS/NS binaries is increased from 346 Mpc to 
364 Mpc, and the NS/NS event rate is increased by a 
factor (364/346)3 = 1.16; see TableU 



assumptions, is 200 Mpc, and the network NS/NS range 
is 200^3 = 346 Mpc. 

Since the only noise source we change, in going from 
one candidate interferometer design to another, is the 
thermoelastic noise, the NS/NS range must be some func- 

tion of [^h(/)/^r(/)]TE = //^BL. 

Buonanno and Chen (private communication) have 
performed the optimization of the advanced LIGO 
optical-noise parameters (the homodyne detection phase 
and the position and reflectivity of the signal recycling 
mirror), as a function of the thermoelastic- noise level, 
to produce for us a curve of optimized NS /NS signal-to- 
noise ratio S/N as a function of {Sh/ S^^)te,- From that 
S / N[{Sh / S^^)te\^ we have computed the corresponding 
network range, (346 Mpc) x {S/N){S/N)^l as a function 
of thermoelastic noise. We show that range in Fig. 0] 

A third measure of a candidate interferometer's per- 
formance is the ratio of its network event rate for NS/NS 
binaries to that of the baseline advanced LIGO network. 
Since the NS/NS binaries are very extragalactic, their 
event rate scales as the range cubed, 

Rate/RatCBL = (Range/346 Mpc)^ . (3.4) 

In a companion paper (15^, OSV evaluate the thermoe- 
lastic noise integral I a numerically for a variety of test- 
mass shapes and beam radii. Here we summarize the 
most important conclusions of those computations and 
their implications for our three performance parameters: 
[Sh{f)/Sl^{f)]T^ [Eq. (El], NS/NS range [Fig.|3|, and 
NS/NS event rate [Eq. (|^ ]. 



3. Cylindrical test masses with MH mirrors and mesa 
beams 

OSV have computed the thermoelastic noise integral 
/ for cylindrical test masses with mesa beams. The test 
masses' volumes were held fixed at the baseline value of 
10** cm^^ (masses fixed at 40 kg), while their physical radii 
Rp and thicknesses H were varied. For each choice of Rp, 
two coated-mirror radii were chosen, R = Rp — 8 mm (the 
baseline choice) and R = Rp (our proposed expansion of 
the coating). In all cases the mesa beam radius D was 
that value for which the diffraction losses are 10 ppm per 
bounce inside the cavity. 

To within the accuracy of their computations, 0.5 
per cent, OSV (Table II of 0) found that the ther- 
moelastic noise integral / is minimized when the test- 
mass dimensions have their baseline values, R = 15.7 
cm, H = 13 cm. In other words, the optimal test-mass 
shape is the same for mesa beams as for Gaussian beams. 
The optimized (10 ppm diffraction loss) radii D for the 
mesa beams, and the values of our three measures of in- 
terferometer performance are shown in Table ^ iii two 
cases: for mirrors coated out to R = Rp — 8 mm (the 
baseline choice), and coated out to R — Rp. 

As is shown in the table, switching from Gaussian 
beams to mesa beams reduces the thermoelastic noise 
Sh oc I by about a factor 3; it increases the NS/NS range 
from 346 Mpc to 465 Mpc if R = Rp - 8 mm, and 497 
Mpc if R — Rp; and it increases the NS/NS event rate 
by a factor (465/346)^ = 2.42 if R = Rp - 8 mm, and to 
(497/364)3 = 2.55 ifR = Rp. 
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TABLE L Optimized light-beam configurations, their ther- 
moelastic noise compared to the basehne (from Table IV of 
fist 'l. their neutron-star binary range, and their event rate 
for NS/NS inspiral divided by the baseline rate. All test 
masses are assumed to be cylinders with the baseline advanced 
LIGO dimensions: physical radius Rp = 15.7 cm and thick- 
ness H — 13.0 cm. The beam radii ro and D are chosen so 
that the diffraction loss per bounce in the arm cavities is 10 
ppm. 



Coated Radius 


Beam Shape 
and Radius 


V ^h^ J TE 


NS/NS 
Range 


Rate 


R — Rp — 8mm 


BL: Gaussian 










To ^ 4.23cm 


1.000 


346 Mpc 


1.00 


R — Rp — 8mm 


mesa 










D/b ^ 3.43 


0.364 


465 Mpc 


2.42 


R^ Rp 


Gaussian 










To — 4.49cm 


0.856 


364 Mpc 


1.16 


R^ Rp 


mesa 










D/b - 3.73 


0.290 


497 Mpc 


2.97 



4- Conical test masses 

By switching from cylindrical test masses to frustums 
of cones, with the same test-mass volume, one could fur- 
ther reduce, substantially, the thermoelastic noise and 
increase the NS/NS range and rate. For detailed explo- 
rations of this, see Tables III and IV and Appendix I of 
OSV d^]. 

We do not discuss this possibility in the present pa- 
per because the current technology for growing sap- 
phire boules, from which to cut the advanced LIGO test 
masses, places a tight limit on the test-mass physical ra- 
dius Rp. It cannot be much larger than the baseline 
Rp — 15.7 cm; and for that maximum radius, and test- 
mass volumes of order the baseline 10** cm^, the optimal 
test-mass shape is cylindrical, with the baseline dimen- 
sions |l5l |. 

When it becomes possible, in the future, to grow 
larger sapphire boules, it might be worth considering test 
masses with frustum shapes p^. 



IV. SENSITIVITY TO MIRROR TILTS, 
DISPLACEMENTS AND FIGURE ERRORS 

The MH mirror figure (Fig. |3J) is somewhat flatter than 
the baseline spherical mirror in its central 10 cm of ra- 
dius where 95 per cent of the light power resides, but 
much more curved in its outer ~ 6 cm. One might 
worry that the central flatness will cause a mesa-beam 
interferometer to be unacceptably sensitive to mirror-tilt- 
induced, mirror-displacement-induced and figure-error- 
induced mixing of parasitic modes into the light beam's 
fundamental, mesa mode. We have investigated this 
mode mixing and find that it is a modest problem, not 
a severe one. We describe our investigations and con- 
clusions in this section. They have been described pre- 
viously in our internal LIGO report |l4j . and a short 



summary of results was given in Sec. II CI above. 



A. Foundations for Investigation 

1. Our tools of analysis 

Our analysis of mode mixing and its consequences is 
based on three independent sets of tools. The first two 
sets are designed for studying the effects of mirror errors 
on the interferometer's high-finesse arm cavities. The 
third set is for analyzing the highly degenerate power- 
recycling and signal-recycling cavities. 

Our first tool set (developed by Richard 
O'Shaughnessy with confirming calculations by Sergey 
Strigin and Sergey Vyatchanin, and described in detail 
in OSV |15| |) is an integral eigenequation for the modes 
of an arm cavity. In the limit of infinite mirror radii (i.e., 
neglecting diffraction losses), the cavity's eigenmodes 
are orthonormal when integrated over the transverse 
plane; this is true for MH mirrors, just as for spherical 
mirrors (Appendix D of [13). OSV have used their 
integral eigenequation to compute the modes with 
untilted, undisplaced and undeformed mirrors and with 
both infinite and finite radii. O'Shaughnessy has then 
tilted, displaced and deformed the ETM of one arm 
cavity and applied first- and second-order perturbation 
theory to its eigenequation to determine the tilt-induced 
and deformation-induced mode-mixing, the resulting 
fundamental eigenmode of the cavity with tilted and 
deformed ETM, its response to the driving beam, and 
the tilt-induced and deformation-induced power going 
out the interferometer's dark port. The details of these 
calculations are given in OSV [l^. In the following 
subsections we describe the main results, we compare 
with computations via our second tool set, and we 
discuss the implications for the use of MH mirrors in 
advanced LIGO interferometers. 

Our second tool set (developed by D'Ambrosio and de- 
scribed in her companion paper [la|) is an adaptation 
of a standard VIRGO and LIGO FFT code for sim- 
ulating interferometers^ (a code originally designed by 
Patrice Hello and Jean- Yves Vinet, then further develped 
by Brett Bochner and others; see and references 
therein). D'Ambrosio has used her adaptation of this 
code to study the same arm-cavity phenomena that OSV 
have studied via the cavity eigenequation and perturba- 
tion theory. She presents the details of her computations 
and some associated perturbation theory analyses in Ref. 
[lfil |. In the following subsections we describe her main 
results, we compare them with the OSV eigenequation 
results, and we discuss their implications. 



^ See http://www.phys.ufl.edu/LIGO/LIGO/STATC.html where 
it is referred to as the FFT Full Field Relaxation Code. 
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Our third tool set (developed by Thorne and described 
in Secs. llV H and HVH below) is designed to deal with the 
influence of mirror errors on the interferometer's power- 
recycling and signal- recycling cavities. Because these 
cavities are highly degenerate, the light trapped in them 
and traversing them is moderately well described by the 
tools of geometric optics, augmented by approximate es- 
timates of the impact of diffraction; so these are the foun- 
dation for our third tool set. 



2. Fiducial mesa and Gaussian configurations 



are those due to diffraction, which we treat as a separate 
mathematical entity from the reflection and transmission 
coefficients. Thus, our power reflection and transmission 
for the ITMs and ETMs are 

rf = 0.995 , tj = 0.005 ; 

r| = 1.0, t| = 0. (4.1) 

We have also assumed the baseline advanced LIGO trans- 
missivities and reflectivities for the power-recycling (PR) 
and signal-recycling (SR) mirrors; see Sees . II V ll and II V Jl 
below. 



When analyzing practical issues, we have carried out 
most of our computations for fiducial MH mirrors that 
have slightly larger coated radii, R]\i — 16cm, than the 
baseline spherical mirrors (14.9 cm at present and 15.7 
cm if our recommendation to coat the outer 8 mm is 
followed), and we have used corresponding /iduciaZ mesa 
beams with modestly larger diffraction losses, Cq — 18 
ppm, than the current baseline of about 10 ppm. In eval- 
uating practical issues in this paper and its companions 
[isL lla | we will compare with fiducial spherical mirrors 
and fiducial Gaussian beams that have this same enlarged 
coated radius Rm = 16 cm and diffraction losses Cq = 18 
ppm. 

Our fiducial mirror and beam parameters, then, are as 
follows: 

• For the fiducial spherical mirrors: mirror radius 
Rm = 16 cm (vs 15.7 cm for the current base- 
line) and Gaussian beam radius ro = 4.70 cm (vs. 
4.23 cm for the current baseline), corresponding to 
a mirror radius of curvature TZc = 83 km (vs. 53.7 
km for the current baseline) and an arm-cavity g- 
factor g = 1 - L/Uc = 0.952 (vs. 0.926 for the 
current baseline). 

• For the MH mirrors: mirror radius Rm = 16 cm 
and mesa beam radius parameter D — Ab — 10.4 
cm, where b = ^/XL/2n = 2.60 cm, with A = 1.064 
/im and L = 4 km the light wavelength and arm 
length. 

[Note: if we had compared the baseline Gaussian-beam 
configuration with a mesa-beam configuration with the 
same diffraction losses, we would have reached approxi- 
mately the same conclusions as we obtain from compar- 
ing these fiducial configurations; see, e.g., the paragraph 
following Eq. (|4.5(l .] 



3. Reflection and transmission coefficients 

For both configurations, fiducial MH and fiducial 
spherical, we have assumed for the ITM mirror the ad- 
vanced LIGO baseline power transmissivity tj = 0.995; 
we have chosen an idealized, perfectly reflecting ETM; 
and for simplicity we have assumed that the only losses 



Driving a Mesa-Beam Interferometer with a 
Gaussian Beam 



One way to produce the desired mesa beams in the 
arm cavities is to drive the interferometer with Gaussian- 
beam light and let the arm cavities or a mode-cleaning 
cavity convert the light into mesa form. One of us, ED 'A 
(Sees. II A and II B of HI), has identified the Gaussian 
beam that has the greatest overlap with the mesa beam of 
an MH-mirrored arm cavity. If one were to drive the MH 
arm cavity directly with a Gaussian beam, this would be 
the driving beam one would want to use. It has a beam 
radius rod = 6.92 cm, compared to our fiducial cavity's 
beam radius at the ITM of = 4.70 cm. This Gaussian 
driving beam Ud{r) has an overlap 



|(uo|u<i)P 



UQUddArea 



0.940 



(4.2) 



with the cavity's fundamental mesa mode uq, which 
means that 94.0 per cent of the Gaussian driving-beam 
light will enter the MH-mirrored cavity, and 6.0 per cent 
will get rejected. See Sees. II A and II B of 16], where 7g 
is denoted C^. 



C. Parasitic Modes in Arm Cavities 

It is useful to think of the MH mirrors as having two 
regions: a central region with radius ~ 10 cm, and an 
outer region extending from radius ~ 10 cm to radius 16 
cm. In its central region, the fiducial MH mirror is much 
flatter than the fiducial spherical mirror; in its outer re- 
gion, it is much more sharply curved; see Fig. |31 The 
flatness of the central region has led to concerns about 
degeneracies of modes and sensitivity to mirror tilts, dis- 
placements and figure errors. 

O'Shaughnessy, Strigin and Vyatchanin have all inde- 
pendently solved the integral eigenequation for the modes 
of a LIGO arm cavity with flducial MH mirrors. They 
have found (Sec. VII A and Table V of T^) that among 
modes that are not strongly damped by diffraction losses, 
the one closest in frequency to the fundamental TEMOO 
mode uq is the lowest TEMOl mode (denoted ui below). 
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Its frequency separation from the fundamental is O.O4O4 
of the cavity's free spectral range, which is 2.5 times 
smaller than for the fiducial spherical-mirrored cavity, 
0.099 X (free spectral range). Evidently, the sharp cur- 
vature of the MH mirrors' outer region compensates suf- 
ficiently for the flatness of their central region, to prevent 
the parasitic modes' frequencies from becoming near- 
degenerate with the fundamental. 



D. Mirror Tilt in Arm Cavities 

Our modeling predicts that mode mixing in the arm 
cavities of a mesa-beam interferometer produced by tilt of 
the ETM's or ITM's should be of no serious consequence, 
if the tilt angles are kept below about 10~^ rad. In the 
following subsections we summarize the calculations that 
lead to this conclusion. 
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1. Parasitic mode mixing in arm cavities 

Two of us have computed the influence of a tilt of the 
ETM on the fundamental mode of an arm cavity: ED'A 
has done this using her EFT code, and RO'S has done 
it by applying perturbation theory to the arm cavity's 
integral eigenequation. The two computations agree on 
the following predictions: 

When the ETM is tilted through an angle 6, the cav- 
ity's fundamental mode gets changed from UQ{r) to 

Uo(r) = (1 - a\/2)un{r) + aiUi(f) + a2U2{r) . (4.3) 

Here r is position in the transverse plane, u„ are unit- 
norm superpositions of modes of the perfectly aligned 
cavity (/ |M„pc?Area = 1), q;„ are mode-mixing coeffi- 
cients that scale as 0", and our computations have been 
carried out only up through quadratic order. The max- 
imum tilt that can be allowed is of order 10~^ radian, 
so we shall express our predictions for the Q!„ in units of 
6*8 = e/lO-^ 

Eor our flducial spherical-mirrored cavities, uq is the 
(0,0) Hermite-Gaussian mode, ui is the (0,1) Hermite- 
Gaussian mode, U2 is the (0,2) mode, and the dominant 
mixing coefficient ai is 

(Eqs. (62) and Appendix F of OSV [H; cf. also Ap- 
pendix A of ED'A 16]). Here g — 0.952 is our fiducial 
arm cavity's g-factor and b = \J XLjl-K — 2.603 cm is its 
transverse diffraction scale. 

Eor the MH cavity, uq is the mesa mode shown (for 
slightly different beam radius) in Eig. [2 and maps of 
Ml, and U2, as computed by O'Shaughnessy [la| by 
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FIG. 5: Contour-diagram maps of the squared moduli |Mi(r)p 
(top) and |u2(j^P (bottom) of the parasitic modes excited 
by tilt of the ETM of a MH arm cavity. These maps were 
generated by O'Shaughnessy fl5| using the eigenequation for 
a cavity with infinite-radius MH mirrors, together with first- 
and second-order perturbation theory. The axis about which 
the ETM is tilted (rotated) is chosen to be vertical (the y 
axis). The units in which |itj(r)|^ is measured are = 
(2.603cm)"^; the peak values are 0.0526 for and 0.0726 

for the peak of |'U2p on the i-axis and 0.0398 for the peak on 
the j/-axis. The contours are separated by 0.01 in units of h~^\ 
the outermost contour in all cases is 0.016~^. See also Figs 
15 and 16 of 0| for maps of these same modes as computed 
by ED'A with her FFT code. 



solving the cavity eigenequation, are shown in Fig. |S1 
D'Ambrosio has used her FFT code to compute maps 
that agree well with these (Figs. 15 and 16 of 16]). The 
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MH mixing coefficients are^ 

af " = 0.02276*8 , a^" = 0.000196i| (4.5) 

(Eqs. (63) of OSV 15]; first equation on page 8 of ED'A 
[l^). Comparison of Eqs. (j4.4|l and H4.5(l shows that, to 
keep the dominant, dipolar mode-mixing coefficient ol\ at 
the same level in the fiducial MH cavity as in the fidu- 
cial spherical- mirrored cavity requires controlling the MH 
mirrors ' tilt angle 3.5 ~ 4 times more tightly. 

[Note: If, instead of comparing our fiducial arm cavities 
(Sec. ITV A 2II . we had compared the baseline spherical- 
mirrored cavity with a MH-mirrored cavity that has the 
same diffraction losses, we would have obtained a^^^ = 
0.0047618 (instead of 0.0064:9^), and^ af^ = 0.0136*8 (in- 
stead of 0.02276*8). We would thereby have concluded 
that, to keep the dipolar mode-mixing at the same level in 
the MH cavity as in the spherical-mirrored cavity requires 
controlhng the MH mirrors' tilt angle 0.013/0.0047 = 
2.8 ~ 3 times more tightly than the spherical mirror's 
tilt (rather than 0.0227/0.0064 3.5 ~ 4 times more 
tightly). This illustrates the fact that the comparison of 
our fiducial configurations gives approximately the same 
answers as would a comparison that focuses on the base- 
line advanccd-LIGO configuration.] 

For our fiducial configurations, the fractional power in 
the dominant, dipolar parasitic mode is of, which will be 
doubled to 

^sph arm total ^ 2«h)2 ^ 0.001(08/3.5)^ , 
pMH arm total ^ 2(af «)2 ^ O.OOl^^ , (4-6) 

when the ITM and ETM are both tilted but about un- 
correlated axes. This suggests that, so far as the arm 
cavities are concerned (i.e., ignoring issues of tilt in the 
recycling cavities), the tilt of MH mirrors need not be 
controlled much better than ~ 10^^ rad. We have ver- 
ified this by examining the effects of the mode mixing 
on various cavity and interferometer parameters; see the 
next three subsections. 



2. Diffraction Losses 

One of us (ED'A), from her EFT simulations, has es- 
timated the influence of ETM tilt on diffraction losses to 
be (last equation on page 9 of [l6| ]) 

C'o"^^ (1 + 0.00491) , (4.7) 

where = 18 ppm for our fiducial mirrors. This result 
has been confirmed to a factor 2 by RO'S by combining 



the clipping approximation with perturbation theory of 
the cavity's eigenequation; Eq. (65) of il^]. The influence 
of ITM tilt should be about the same, thus doubling the 
coefficient of This tilt-induced increase in losses is so 
small that it can be ignored for tilt angles below ~ 10~® 
rad. 



3. Arm Cavity Gain 

We have computed (ED'A via the EET code and RO'S 
via perturbation theory) the following inffuence of ETM 
tilt on the arm cavity gain*" 

Gf^ = 740(1 - 0.000576*8) (4-8) 

(Eq. (68) of [13; Sec. HI B of 13). This result assumes 
the baseline values for the power transmissivities of the 
ITM and ETM, and assumes for simplicity that the only 
losses are diffraction losses; Sec. IIV A 31 The factor 740 
assumes the cavity is driven by its best-fit Gaussian beam 
and thus is smaller by about 7g = 0.940 than the gain in 
the untilted, fiducial spherical-mirrored cavity. The tilt 
of the ITM should produce about the same gain reduction 
as that of the ETM, thus doubling the coefficient of 6*| 
to ~ 0.001. This coefficient is small enough that the tilt- 
induced decrease of MH arm-cavity gain will be negligible 
if 9 is controlled to ~ 10^^ rad. 



4- Dark Port Power 

We have computed the influence of the tilt of one ETM 
on the dark-port output light (ED'A using the EFT code 
and RO'S using perturbation theory). Multiplying that 
result by four to deal with the case of all four cavity mir- 
rors being tilted about uncorrelated axes, we find for the 
fraction of the interferometer's input power that winds 
up at the dark port in the fundamental mode uq and the 
parasitic modes ui and U2'. 

pMH DP total ^ ^ Qgi ^ 

pMH DP total ^ 472(af H)2 ^ 2OOO0| ppm , 

pMH DP total ^ Aj^iaf'f ~ O.Wj ppm . (4.9) 

(Eqs. (70) of ^3 multiplied by 4; last four equations on 
page 8 of |16|^ multiplied by 2). Without an output mode 
cleaner, the dark-port power would primarily be in the 
dipolar mode ui, and for 9 < 10^^ rad it would constitute 
< 0.2 per cent of the input light. The planned output 
mode cleaner will wipe out this u\ power and the power in 



5 RO'S and ED'A both find af^ = 0.022798; for af^ RO'S finds 

0.000186»| and ED'A finds 0.00020e|. 
^ We thanlc Pavlin Savov for computing this number for us. 



RO'S finds = 737(1 - 0.000556»|); ED'A finds ^ 

746(1 - O.OOO590|). 
8 In Hi Pdp/Po is half our P^H DP total ^ P^°"'i'P is half 
our (Pg^^ -|- P|^^ totalj 
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mode U2, leaving only the tiny fundamental-mode power, 
which should be totally negligible for below 10~* rad. 

For comparison, the dark-port powers with the fiducial 
Gaussian beams are 

psph DP total ^ 4(c,«Ph)4 ^ 1.0(08/3.5)4 ppm , 
psph DP total ^ 4(^-^Ph)2 _ 2000(^?8/3.5)2 ppm , 

(4.10) 

which shows once again that the fiducial spherical- 
mirrored arm cavities are ~ 4 times less sensitive to tilt 
than the fiducial MH arm cavities. 



E. Transverse Displacement of Arm Cavities' 
Mirrors 

When the ETM is displaced transversely through a 
distance s, the cavity's fundamental mode gets changed 
from Uo{r) to 

<W = (1 - Ci/^)Mr) + Ci^iW + C2W'2(r) , (4.11) 

where the parasitic modes i(j„, like the m„'s, have unit 
norm, {wn\wn) = J |w„pdArea = 1, and have phase ad- 
justed so the coupling coefficients C„ are real, and where 

R'OS has computed the coupHng coefficients Cn for 
fiducial spherical mirrors and fiducial MH mirrors by ap- 
plying perturbation theory to the cavity's eigenequation; 
and for MH mirrors, ED'A (unpublished) has computed 
the influence of the ETM displacement on the arm cav- 
ity's gain using her FFT code, and has then deduced Ci*^ 
by fitting to an analog of Eq. (68) of J^J in the limit of 
small displacements. For both spherical and MH mirrors, 
^2 is negligible compared to Ci when the displacement is 
s ^ 6 = XL/2tt = 2.60 cm, so we shall ignore C2- 
From the R'OS and ED'A computations, we deduce that 
the fiducial MH-mirrored cavities and fiducial spherical- 
mirrored cavities are approximately equally sensitive to 
transverse displacements of the ETM; their coupling co- 
efficients are:^ 

= 0.010s„„„. (4.12) 

Here g = 0.952 is the fiducial arm cavity's g-factor, and 
Smm is the ETM's transverse displacement in millime- 
ters. The MH cavity's lesser sensitivity to displacement 



9 Here Ci^*" and (f^^ are from R'OS 15| Eqs. (57) and (58). ED'A 
(unpublished) obtains Cj^H ~ 0.008s mm ; 20 per cent smaller than 
R'OS; we regard her Ci ^ as less reliable than R'OS's because 
the phase of a parameter 72 used in her fit is not known with 
sufficient reliability. 



presumably arises from the flatness of its mirrors in the 
central region, where most of the light power resides. 

The corresponding fraction of the arm-cavity carrier 
power driven into the (dipolar) parasitic field wi is (cf. 
Eq. (59) of 05.]) 

parm total = /•2 ^ f 100(s/1.3mm)2ppm spherical, 
^ I 100(s/1.0mm)2ppm MH. 

(4.13) 

The fraction of the input carrier power driven out the 
dark port when the ETMs of both arm cavities are dis- 
placed through a distance s but in uncorrelated directions 
is about twice the above [cf. J^], Eq. (60) multiplied by 
two, since here both ETMs are displaced and there just 
one]: 

pDP total ^ OtV^ ~ / 190(s/1.3 mm)2ppm spherical, 
1 7o<,i - ^ 190(5/1.0 mm)2ppm MH. 

(4.14) 

These coupling coefficients and parasitic-mode powers 
are sufficiently small that transverse displacements are 
not a serious issue, and so shall ignore them in the rest 
of this paper. In any event, the low sensitivity to a change 
from fiducial spherical to fidcuial MH mirrors makes dis- 
placement a non-issue in the any decision about whether 
to use MH mirrors. 



F. Errors in the Arm Cavities' Mirror Figures 

1. Billingsley's Worst-Case Figure Error 

Garilynn Billingsley (of the LIGO Laboratory, Cal- 
tech) has provided us with a map of a worst-case figure 
error, Szyjci^TV) [height error as function of Cartesian 
coordinates in the transverse plane], produced by cur- 
rent technologies. Her map is based on the measured 
deviation of a LIGO-I beam-splitter substrate from flat- 
ness. The measured substrate had diameter 25 cm; she 
stretched its deviation from flatness (its "figure map" ) to 
the fiducial mirror diameter of 32 cm, fit Zcrnike poly- 
nomials to the stretched map, and smoothed the map by 
keeping only the lowest 36 Zernikes. 

We show a contour diagram of the resulting figure map 
(figure "error" ) in Fig.|H| In the central region (innermost 
10 cm in radius), the peak to valley error Az is about 30 
nm, while in the outer region (10 cm to 16 cm in radius), 
it is about 110 nm. Billingsley thinks it likely that in 
the central region (which dominates our considerations), 
peak-to- valley errors of Az ~ 5 nm may be achievable for 
MH mirrors — about 1/5 as large as in Fig. El and we 
have found that Az — 0.2 x 30 = 6 nm is small enough 
that the influences of the figure error scale, for Az < 6 
nm, as Az or Az^ with higher-order terms producing 
< 10 per cent corrections. Accordingly, in the analyses 
described below we shall use Billingsley's map, scaled 
down in height by a factor e: 

Sz = eSz„cix,y) , (4.15) 
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FIG. 6: Contour diagram of Billingsley's worst-case figure 
error [height Szv,c in nanometers as a function of transverse 
Cartesian coordinates {x,y) in centimeters]. The hash at the 
outer edge of the mirror is an artifact of our numerical ma- 
nipulation of Bilhngsley's map. 

and we shall use e = 0.2 and Az = 6 nm as our fiducial 
values for e and Az. Jean Marie Mackowski (an expert 
in coating mirrors) believes that Az ^ 2 nm errors may 
be achievable for MH mirrors, if the mirror figure is pro- 
duced by coating; this corresponds to £ ~ 0.07. For 
comparison, the figure errors for the spherical LIGO-I 
mirrors are Az ~ 3 to 6 nm (rms errors ^ 1 to 2 nm). 

2. Mode Mixing by Figure Errors Without Compensating 

Tilt 

We have computed the mode mixing in an MH cavity 
when Billingsley's worst-case figure error, reduced by e ^ 
0.2, is placed on the ETM. As in our tilt studies, the 
computation was done independently by ED 'A using the 
FFT code and by RO'S using arm-cavity perturbation 
theory. 

By analogy with Eq. (|4.3|) , the fundamental mode with 
deformed ETM has the following form 

u'o = il-f3f/2)uo + l3ivi, (4.16) 

where the parasitic mode wi, like the u„'s, has unit norm 
{vi\vi) — J juipdArea — 1, and has its phase adjusted so 
that Pi (X e/0.2 is real. By contrast with the tilt-induced 
mode mixing, where ui is dipolar (angular dependence 
cos(^), the deformation parasite vi has a complicated 
shape that depends on the details of the deformation and 
that therefore contains a number of multipoles. A map 
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FIG. 7: The power distribution |ito — ?ioP = (in units 

1/m^) of the deformation-induced parasitic mode when the 
deformation e5zwc(a;,2/) with e — 0.2 is applied to the ETM 
of an MH arm cavity. This map was computed by ED'A 
1^ using the FFT code; the map computed by RO'S using 
perturbation theory fl^ is in reasonable agreement with this 
one (e.g., the heights of the two peaks are {0.111,0.080} in 
the FFT map, and {0.110,0.088} in the perturbation map). 
The parasitic power |/9iwi|^ scales as (e/0.2)^. 

of the power distribution of the admixed mode is 

shown in Fig. \7\ 

The fraction of the arm cavity power in the parasitic 
mode is 

pMH arm ^ ^2 ^ o.0011(£/0.2)2 (4.17) 

(Eqs. (73) of ^3; first term in first equation on page 
10, and also first equation on page 11, of 0|). The 
fraction of the interferometer's input power that goes out 
the dark port (if the interferometer is driven by the best- 
fit Gaussian mode and if only one of the arm mirrors 
— one ETM — is deformed) is given by 

pMH DP ^ ^2^2 ^ o.0010(£/0.2)2 (4.18) 

(Eq. (75) of O; cf. also O). At the leading, e^, order 
in the deformation (the order to which we have com- 
puted), this dark-port power is entirely in the parasitic 
mode vi . Our FFT and perturbation-theory calculations 
agree on the parasitic powers H4.17|l and (|4.18|) to within 
about five per cent. 

The dark-port power H4.18|l and parasitic arm-cavity 
power (|4.17() are influenced primarily by the figure error 
in the central (10 cm radius) region of the ETM, because 
about 96 per cent of the mesa-mode power is contained 
in that central region, and only about 4 per cent in the 
outer region — and of the outer 4 per cent, 3/4 (3 per 
cent) is in the annulus between 10 and 11 cm. The insen- 
sitivity to outer-region deformations is fortunate, because 
Billingsley tells us that it will be much easier to keep the 
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figure errors small in the central region than in the outer 
region. 

[We have verified the insensitivity to the outer-region 
deformations by evaluating (R'OS via perturbation the- 
ory and ED 'A via FFT simulations) the dark-port power 
for a mirror deformation of the form 



5z = - 



eo^r". (4.19) 

Here Jz™"'"''*' is equal to 5z^c at r < 9.6 cm and is zero at 
r > 12.2 cm, and between 9.6 and 12.2 cm, Sz^^^""^ / Sz^^ 
falls linearly from 1 to zero; and similarly Sz'^^'^ is equal 
to J^wc at r > 12.2 cm and is zero at r < 9.6 cm, and 
between 9.6 and 12.2 cm 5z°^*®''/5zwc grows linearly from 
to 1. We find, as a function of the central-region and 
outer-region weightings, ^° 

pMH DP ^ ^2^2 ^ o.0010[(e,/0.2)2 + [eo/O.lf] ; 

(4.20) 

so that, for example, if Billingsley's worst-case perturba- 
tions are reduced by tc = 0.2 in the central region (to 
Az — 6 nm), but are are reduced only to to — 0.7 in the 
outer region (so Az = 21 nm there), the outer region will 
contribute about as much power to the dark port as the 
inner region.] 

When all four arm-cavity mirrors are subjected to un- 
correlated deformations, the arm-cavity parasitic power 
(|4.17l) will be increased by a factor 2 and the dark-port 
power (|4.18|l by a factor 4, to 

pMH arm total ^ ~ 0.0025( Az/6nm)2 , 
pMH DP total ^ 4^2^2 _ 0.005 (Az/6nm)2 (4.21) 

where Az is the peak-to- valley mirror deformation in the 
central region. This suggests that, so far as arm-cavity 
mode mixing is concerned, peak-to- valley figure errors of 
order 6 nm in the inner 10 cm are acceptable. 



3. Mode Mixing by Figure Errors With Compensating Tilt 

The parasitic mode vi (Fig. [7J contains a significant 
amount of dipolar field, as one sees from the asymme- 
try of the map. The advanced LIGO tilt control system, 



10 The second coefficient, (eo/0.7)2, is from the ED'A FFT calcu- 
lation s: a fit to numbers in the first and third equations on page 
If of If6l . in which fit we (arbitrarily) ignored the possibility 
of a cross term between to and tc- The RO'S perturbation- 
theory calcuations use mirrors of infinite radius, and the modes 
admixed by the outer-region perturbations have sufficient am- 
plitude outside the radius of our fiducial mirror, r > 16 cm, 
as to make his results for the dependence on to somewhat un- 
reliable. (The unrelieability arises from two effects: his ad- 
mixed modes are wrong due to omitting the influence of the 
mirror edge, and his analysis misses the influence of the mir- 
ror errors on diffraction losses.) His analog of Eq. 14.201 is 
pMHDP ^ 0.0010[(ec/0.2)2 + e„(e„ - tc)/{0.7)% implying 
5 times weaker sensitivity to outer-region perturbations than 
central-region, by contrast with the 3.5 times weaker sensitiv- 
ity in ED'A's more reliable Eq. l4.2Ui . 
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FIG. 8: Billingsley's worst-case figure error reduced by a 
factor e = 0.2, when the mirror is tilted to optimally re- 
duce the mesa beam's odd-parity power: 5z — 0.25zwc -I- 
Scomponsatc sin((/3 — 55°)]. The numbers on the contours are 
height in nanometers. 



based on a quadrant-diode readout of asymmetry in the 
power distribution u'q, will tilt the mirror so as to remove 
the overlap between the deformed parasitic field vi and 
the dipolar-tilt parasitic field ui. ED'A and R'OS have 
independently computed that the optimal tilt is about 
^'compensate = 1-3 X 10~®(e/0.2) radiaus about a line ro- 
tated 55 degrees from the x axis (Eqs. (77) of O), and 
have computed the resulting field u'q = uq + PiVi + aiui 
with minimum parasitic-mode power. Figure|Slshows the 
mirror deformation after tilt, and Fig.|^shows the para- 
sitic power distribution |/3iWi +q;iUi p for e — 0.2. Notice 
that the tilt has largely but not completely removed the 
dipolar asymmetry. Some residual dipolar field remains 
— that portion whose radial distribution prevents it from 
being compensated by a tilt. 

We denote by vic the tilt-compensated parasitic mode 
and by /?ic its (real) amplitude, so PicVic = PiVi +aiui. 
Then the cavity's eigenmode, with tilt compensation (in- 
cluding the second-order loss of power from uq) is 



u'q = {1 - Pfj2)uo + PlcVlc 



(4.22) 



and optimization of the tilt compensation gives for the 
fraction of the power in the parasitic mode 



pMH arm ^ ^2^ ^ 0.00040(e/0.2)2 = 0.00040(Az/6nm)2 

(4.23) 

(Eq. (79) of 15]; last equation on page lOof [ig, divided 
by 7q — 0.940). The fraction of the interferometer's input 
power that goes out the dark port (all in the parasitic 
modes), after this tilt compensation, is 

pMH DP ^ ^2^2^ ^ 0.00038(e/0.2)2 (4.24) 
(Eq. (80) of l^st equation on page 10 of 16]). 
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FIG. 9: The power distribution \u'q — uo\'^ — \/3icVic\'^ (in units 
1/m^) of the deformation-induced parasitic mode when the 
tilt-adjusted deformation Sz = 0.2(5 Zwc + ^compensate sm{if — 
55°)] is applied to the ETM of an MH arm cavity. This map 
was computed by RO'S by applying perturbation theory to 
the cavity's eigenmode llSll . The map computed by ED'A 
using the FFT code fl^ agrees reasonably well. 



Comparing with Eq. (|4.18ll . we see that the figure-error- 
induced power out the dark port is reduced by a factor 3 
by the compensating tilt of the deformed mirror. 

The dark-port power (|4.24() is for an interferometer 
with one deformed MH arm mirror. When all four mir- 
rors are deformed in uncorrelated ways, the arm-cavity 
parasitic power will be doubled and the dark-port power 
will be quadrupled: 



MH arm total 



pMH DP total 
^1 



2Pl^ ~ 0.0008(Az/6nm)2 , 
47^/32^ ~ 0.0015(Az/6nm)2 (4.25) 



This suggests that, so far as arm- cavity mode mixing is 
concerned, we could live with central-region deformations 
as large as 10 nm. Recycling-cavity issues, discussed be- 
low, will place much tighter constraints on the mirror 
figures. 

We note in passing that the computations reported 
in this subsection have verified the ability of the LIGO 
control system to measure the lowest-radial-order dipolar 
component of a mirror figure error, and remove it by 
applying a tilt. 



4. Comparison with Mode Mixing for Baseline Spherical 
Mirrors with Tilt-Compensated Figure Errors 

As part of her modeling of the baseline configura- 
tion of advanced LIGO (a different project from that re- 
ported here), Ed' A has evaluated tilt compensation and 



mode mixing for an interferometer with the baseline (BL) 
spherical mirrors, to which the figure error of Fig. |^ with 
£ = 0.2 has been applied. For the BL spherical mirrors, 
95 per cent of the light power is contained inside a circle 
with radius ^95% ~ 7.3 cm, by contrast with the fiducial 
MH mirrors, where ^95% ~ 10 cm (and the fiducial spher- 
ical mirrors with rg5% ~ 8 cm). Since (7.3/10)^ ~ 0.5, 
the BL Gaussian beam feels a central region of the de- 
formed mirror with area approximately half that felt by 
the fiducial mesa beam. 

Interestingly, the tilt compensation that removes the 
lowest-radial-order dipolar deformation is twice as large 
for the BL spherical case as for the fiducial MH case: 
C„;lsatc ^ 2.6x 10-8(e/0.2) compared to ^^4" ^^^^ ^ 
1.3 X 10"^(£/0.2). The BL tilt is about a line rotated 
76 degrees from the x-axis, compared to 55 degrees for 
the fiducial MH case. These differences are due to the 
different regions of the mirror that the two beams feel. 

After tilt compensation, the peak-to- valley figure-error 
variations in the 95-per-cent-power region are [for figure 
errors reduced from worst case by e = 0.2] (Az)gg^^'^'^ ~ 



fid MH 

95% 
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2 nm (inner 7.3 cm), by contrast with (Az) 
nm (inner 10 cm). This factor 3 difference is caused by 
Billingsley's figure errors becoming larger rapidly with 
increasing radius — an issue that will have to be dealt 
with if MH mirrors are used in advanced LIGO. 

D'Ambrosio's analysis of mode mixing by the tilt- 
compensated figure errors gives, for the fraction 
of arm-cavity power driven into parasitic modes, 
pBLsph arm total _ 0.0008 [(Az)95%/12nm] ^ , Compared 



to Pi 



MH 



I total 



G.0008[(Az)95%/6r 



[Eq. 



Thus, for the same amount of mode mixing, the BL 
spherical mirrors can have twice as large peak-to-valley 
figure errors as the fiducial MH mirrors, and in the BL 
spherical case the relevant figure errors are confined to 
a region half as large in area {rg^% ~ 7.3 cm) as in the 
fiducial MH case (^95% ~ 10 cm). Presumably, the two- 
fold greater sensitivity of MH arm cavities to figure errors 
is caused by the greater degeneracy of their modes (last 
paragraph of Sec. lIV CJ|l . If we had compared fiducial MH 
mirrors with fiducial spherical mirrors instead of with BL 
spherical mirrors, the diffferences would have been more 
modest. 



G. Influence of Mirror Tilt and Figure Errors on 
Thermoelastic Noise 



When one MH mirror of an arm cavity is given the 
deformation e^Zwc, the resulting deformation of the arm- 
cavity power distribution, (5|uq|2 = [mqI^ — |uo|^ increases 
the thermoelastic noise. The following argument (due to 
RO 'S Sees. IV C 4 and VII E 4 of O ) shows that, at lead- 
ing (linear) order in e/0.2, only the circularly symmetric 
portion of the parasitic mode (3iVi — u'q — uq contributes 
to the thermoelastic noise increase: The fractional noise 
increase is equal to the fractional increase in the ther- 
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moelastic noise integral I = J (V0)^(ivolume [Eqs. (|3.1() 
and 

rc-TEMH ST r ^ 

= y cx y (Ve) • (V5e)dvolume . (4.26) 

Here Q is not a tilt angle but rather is the expansion 
(fractional increase of volume) of the substrate mate- 
rial when a static pressure P cx |mqP is applied to the 
mirror face. Since the unperturbed mesa beam is circu- 
larly symmetric, so will be the unperturbed expansion 
Q, which means that only the circularly symmetric por- 
tion of the expansion perturbation SQ, and thence the 
circularly symmetric portion of the pressure perturba- 
tion SP (X 5|moP will contribute to the noise increase. At 
leading (linear) order in the mirror deformation e/0. 2, the 
circularly symmetric portion of <^|uoP arises solely from 
the circularly symmetric portion of Su'q = PiVi; thus, as 
claimed, only the circularly symmetric portion of Pivi 
can increase the thermoelastic noise. 

This same argument shows that the mesa-mode de- 
formation produced by mirror tilt cannot influence the 
thermoelastic noise at first order in the tilt angle; and 
therefore, we need not be concerned about the influence 
of tilt on the thermoelastic noise — whether the tilt is 
unintended, or is being used in a controlled way to com- 
pensate the errors in the mirror figures. 

RO'S and SS have independently computed the frac- 
tional increase in the thermoelastic noise integral at lin- 
ear order in e/0.2; cf. Eq. (81) in OSV Their result, 
multiplied by four to account for four arm-cavity mirrors, 
is 

r oTE, MH total 

' TE.MH = 0.14(e/0.2) c 0.14(Az/6nm) . (4.27) 

This 14 per cent increase of S^, when all four mirrors 
are subjected to 6 nm figure errors in their central re- 
gions, is to be compared with the factor 1/0.34 = 295 per 
cent decrease in thermoelastic noise achieved by switch- 
ing from spherical mirrors to MH mirrors. There may 
also be a small increase in thermoelastic noise when 
a spherical mirror is deformed. Assuming, conserva- 
tively, no deformation-induced noise increase for spher- 
ical mirrors, switching from spherical to 6nm-deformed 
MH mirrors will reduce the thermoelastic noise by a fac- 
tor 0.34 X 1.14 = 0.39, which in turn will increase the 
distance for NS/NS binaries from 364 Mpc to about 490 
Mpc (Fig.^ and increase the event rate by about a factor 
2.45; cf. Table D 

H. Mesa vs. Gaussian Beams in Recycling Cavities 



cavities, which strongly enhances error-induced mode 
mixing. One might worry that for MH mirrors, with their 
greater central-region flatness, this severe mode mixing 
might be made substantially worse. 

We have examined this question and conclude that 
for the two wideband advanced LIGO interferometers, 
there is not much difference between the fiducial Gaus- 
sian beams and the fiducial mesa beams, with regard to 
their susceptibility to mode mixing in the recycling cavi- 
ties. The only significant difference arises from the fact 
that the mesa beams are larger and therefore sample, 
with significant power, larger-radii regions of the mirrors 
(the regions between, say, 8 cm radius and 10 cm radius), 
where the deformations may be worse. 

The reason that the mode mixing is only marginally 
sensitive to the beam shape is quite simple: If hetero- 
dyne readout were used, then once RF modulated side- 
band light gets into the RF optical cavity bounded by the 
power recycling (PR) mirror, the signal recycling (SR) 
mirror and the two input test-mass (ITM) mirrors, the 
sideband light would make roughly A/rf ^ ^(cavity fi- 
nesse J^pr) ~ (7r/2)/(l — Rpr) ~ 25 round trips before 
losing 95 per cent of its power out the beam splitter's 
dark port.^^ (Here i?pR = ii!sR — 0.94 is the power 
reflectivity of the PR and SR mirrors; they are made 
equal so as to transmit the RF light with maximum ef- 
ficiency to the photodetector.) And once signal light 
gets into the signal recycling cavity, it makes roughly 
A/sR ~ \^ /{^ ~ pVR) ~ 40 round trips before losing 95 
per cent of its power out the dark port or back into the 
arm cavities. (Here R = 0.995 is the ITM power reflec- 
tivity and p = \/0.93 is the amplitude reflectivity of the 
SR mirror, in the notation of Buonanno and Chen Q.) 
The Fresnel length (transverse diffraction scale) for light 
that makes Af round trips in either recycling cavity, with 
cavity lengths inp — 20 m and isR — 10 m, is 

rp = ^/Xjmf ~ 3cm , (4.28) 

for both cavities; here Ao — l/xm is the light's wavelength. 
This Fresnel length is ~ 1/2 of the ~ 5 cm scales on which 
the ideal mirror shapes and the central-region worst-case 
mirror errors vary, and it is small compared to the ~ 15 
to 20 cm diameter beams themselves. There thus is only 
modest diffractive coupling between light rays, and the 
light bouncing back and forth in each recycling cavity is 
describable, to moderately good accuracy, by geometric 
optics. Moreover, because the mirrors (whether MH or 
spherical) are nearly flat and nearly identical, the light's 
rays, to rather good accuracy, are all parallel to the optic 
axis and to each other and are thus decoupled from each 
other. If the mirrors and beam splitter were perfect and 



In LIGO-I interferometers, and in the baseline design 
for advanced LIGO interferometers, the greatest sensitiv- 
ity to mirror tilts and figure errors occurs in the power re- 
cycling (PR) and signal recycling (SR) cavities. This sen- 
sitivity arises from the near degeneracy of the recycling 



Our definition M = (number of round trips) ~ .7^/2, based on 95 
per cent power loss, differs by 4/7r ~ 1.27 from a more conven- 
tional M = (2/it)J^ based on the phase shift of light emerging 
from the cavity. 
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ideal in shape, the extreme length of their radii of curva- 
ture, > 50 km, compared to the optical pathlength in the 
recycling cavities, 2iM ~ 1 km, would guarantee that the 
mesa beam would resonate equally well in the ideal MH- 
mirrored cavity or in the ideal spherical-mirrored cavity, 
or in a precisely flat-mirrored cavity; and the flducial 
Gaussian beam would also resonate, equally well, in all 
these cavities. 

If an advanced LIGO interferometer is operated in 
narrow-band mode, then the number of round trips the 
signal light makes in the SR cavity will be much larger 
than 40, and the geometric optics approximation will be- 
gin to fail significantly. More specifically, for ITM re- 
flectivity R = 0.995 and optimized narrow banding at 
{500 Hz, 1000 Hz}, the SR mirror's amplitude reflectivi- 
ties are p = {0.994,0.9985} [see discussion following Eq. 
(|4.35|) below] , corresponding to a number of round trips 
in the SR cavity Wsr ~ {180, 400} and Fresnel lengths 
rp ~ {6cm, 9cm}. These Fresnel lengths (the transverse 
scale for diffractive fight spreading) are about 1/3 to 1/2 
the 95-percent-power diameter of the beam, 16 cm (fidu- 
cial Gaussian) and 20 cm (fiducial mesa). As we shall 
see below, this means that geometric optics can be used 
to get a rough upper limit on the fractional decrease of 
signal strength due to tilt and irregularities of the SR 
cavity's mirrors. 

I. Decrease of Signal Strength Due to Mirror Tilts 

1. Foundations 

The mirror tilts produce a mismatch between various 
modes of the light, thereby decreasing the signal strength 
at the interferometer output (the beam splitter). We 
shall focus on the signal-strength decrease at the fre- 
quency of the signal light's optical resonance in the arm 
cavity. This optical resonance is the one that is used 
to produce a noise [Sh {f)] minimum for a narrowbanded 
advanced LIGO interferometer, and it is the right-hand 
minimum (at / ~ 230 Hz) of the optical noise for the 
standard wide-band advanced LIGO interferometer (Fig. 
0, and approximately the minimum of the wide-band 
interferometer's total noise. At this frequency, for MH 
mirrors with their reduced thermoelastic noise, the dom- 
inant noise source is shot noise; for spherical mirrors, 
short noise and thermoelastic noise are comparable at 
this frequency, but in a small piece of our analysis, we 
shall idealize the shot noise as dominant. 

The (unit-norm) modes whose mismatch decreases the 
signal strength are the following: 

1. Mo, the eigenmode of perfect arm cavities. 

2. Uq, the carrier's eigenmode in an arm cavity with 
tilted ITM and ETM. 

3. Ug, the signal field's eigenmode in an arm cavity, at 
the center of its optical resonance, with tilted ITM, 



ETM, and SRM (signal recycling mirror). 

4. U(, the field produced when uq is transmitted 
through the signal recycling cavity with tilted ITM 
and SRM. 

5. u'j., the reference- light field that is beat against the 
signal light to produce the input to the photode- 
tector. For the baseline homodyne readout system, 
this will be carrier light u'q transmitted through the 
ITM and SR mirror to the photodiode. If hetero- 
dyne readout were used, it would be side-band light 
transmitted through the RF cavity to the photode- 
tector. 

For each primed field u'^, we denote by Si the fraction 
of its light power that is in parasitic modes and thus has 
been lost from the fundamental mode uq due to mirror 
tilt: 

{uo,u'i)^ (4.29) 

The signal amplitude entering the photodetector is 
proportional to 

5« «,f'0«,u;,)«,Md) . (4.30) 

The sequence of terms, from right to left, have the fol- 
lowing meanings, and we approximate them as follows: 

1. {u'q, Ud) describes the influence of tilts on the driv- 
ing of the arm cavity's eigenmode by the Gaussian 
driving fleld. For simplicity, we neglect the tiny 
coupling of the driving field to the second-order 
parasitic mode U2 contained in uq and therefore 
approximate this coupling amplitude by {uQ,Ud) — 

7o(m[„mo) =7o(l-afB/2-ay2) =7o(l- V2). 
Here the subscripts I and E denote the contribu- 
tions from the tilts of the ITM and ETM. 

2. (u(,,Mq) describes the influence of tilts on the driv- 
ing of the arm cavity's signal fleld by its car- 
rier field (via the gravitational-wave-induced mo- 
tion of the mirrors). For simplicity we neglect the 
(nonzero) overlap between the parasitic modes con- 
tained in u'^ and u'q, thereby obtaining (m'j,,Mq) — 
(l-J,/2-Jo/2). 

3. {u'j.,f'u'g) describes the influence of tilts on (i) the 
passage of the signal u'^, through the SR cavity 
(with cavity transmissivity f in the notation of 
Buonanno and Ghen '^), and on (ii) the over- 
lap of the transmitted signal light with the refer- 
ence light to produce the photodetector current. 
Again we neglect correlations between the para- 
sitic components of the fields and therefore ap- 
proximate the infiuence of the tilts by (m^, fu'^) oc 
{1 — Sr/2 — Ss/2 — 6t/2). The Sj terms represent the 
loss of overlap due to the parasitic-mode fields (as- 
sumed uncorrelated) contained in u'^ (the Sr term) , 
contained in u'^ (the Sg term), and generated by 
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the passage of the signal light through the SR cav- 
ity, whose mirror tilts deform the transmissivity f', 
(the 6t term). 

If there is no mode cleaner on the interferometer out- 
put, then the rms amplitude of the shot noise (the dom- 
inant noise source at our chosen frequency) is A'^ oc 
^ {u'^, u'^) = 1; i.e. the parasitic- mode components of 
mJ, contribute to the rms noise amplitude along with 
the fundamental- mode component. However, a mode 
cleaner will remove the parasitic components, so that 

iV (X (wcO = 1 - <5r/2. 

Combining the above approximations to the various 
terms, we find for the ratio of noise power to signal power 
(which is proportional to the spectral density of noise 
Sh referred to the gravitational- wave input h), at the 
minimum of the optical resonance: 

S (X (xl + 25 +2S + (5 + I ' mode cleaner 
jS'p ° * \ with mode cleaner. 

(4.31) 

(Note that we have ignored the increase in shot noise due 
to carrier-light parasitic fields going out the dark port, 
Eqs. H4.9|l and (|4.25l) . under the assumption that it is 
negligible, either because of an output mode cleaner or 
because the arm-cavity-mirror figures and tilts are ade- 
quately controlled.) 

We shall now examine the various contributions to the 
increase of 5*^, one by one. 



2. Carrier Light in Arm Cavity 

The fraction of the arm-cavity carrier power that is 
driven into parasitic modes by tilts of the ETM and ITM 
a^j. The loss of this carrier power to par- 



is <^o = 

asites reduces the signal power and thence increases 
the gravitational-wave noise by 



SSh 
Sh 



MH 



= 2do = 4ai = 0.01 



^2 X 10-8^ 

(4.32) 

Here we have assumed that both mirrors are tilted 
through the same angle 6 but about uncorrelated axes, 
we have assumed MH mirrors, and we have used Eq. H4.5|l 
for ai. The fiducial spherical mirrors are four times less 
sensitive to tilt, so to keep this contribution to Sh below 
one per cent, we must control the ITM and ETM tilts to 
an accuracy 



qMH 

^1% 



= 2 X 10" 



91^^ = 8 X 10-« 



(4.33) 



These are modest constraints on tilt. 



3. Signal Light in Arm Cavity 

The signal recycling (SR) cavity presents a complex 
amplitude reflectivity p' = e~'^l^ e^^^^ to the arm cav- 
ity's signal light (Eqs. (5) and (13) of Buonanno and 



Chen 13). Here F = c/2L is the interferometer's free 
spectral range, e — e(i?, p, 0) and A = A(i?, p, </)) are real 
functions of the ITM power reflectivity _R, the SRM am- 
plitude reflectivity p and the SR cavity's tuning phase 
(/i — (fc^)mod27r, with t the length of the cavity; and our 
notation is that of Buonanno and Chen . Tilts of the 
ITM and SRM produce a spatially variable reflectivity 
p' . The spatial variations of the modulus e~'^l^ of p' pre- 
sumably will have much less influence on the arm cav- 
ity's signal eigenmode than the spatial variations of 
the phase. (This claim deserves to be checked.) Assum- 
ing this is so, then the dominant influence of an ITM 
or SRM tilt d is to produce a spatially variable mirror 
displacement 



8z = 0r sin ip 



(4.34) 



(where ip is azimuthal angle and r is radius), which 
in turn (in the SR cavity's geometric optics limit) pro- 
duces a spatially variable phase of the cavity reflectivity, 
arg(p') = SX/F = {dX/d(j)){k/ F)Sz. If the cavity were 
replaced by a single mirror that is displaced through a 
distance Szcs, then this phase change would be 2k6zcs. 
Correspondingly, the tilt of the ITM or SRM produces 
an effective mirror displacement ^Zoff = ASz, where the 
amplification factor A is given by 



A 



5z^, 



d\/d(j) 
2F 



(1 - R)p 



2p+ (1 + p2)cos20 



(1 + p2) +2/9COs2(/) 

(4.35) 

see Eq. (18) of Buonanno and Chen Q- 

We shall focus on three configurations for the SR cav- 
ity: (i) The standard wideband advanced LIGO configu- 
ration (denoted "WB"), for which R = 0.995, p = V0J3, 
and (j) — 7r/2 — 0.06. (ii) An interferometer narrowbanded 
at a frequency / — A/27r ~ 500 Hz with bandwidth 
A/ = e/27r ~ 50 Hz, which has a noise minimum of 
~ 1 X IO^^'^'/VHz; this configuration (which we shall 
denote "500") is produced by i? = 0.995, p = 0.994, 
and (j) — 1.541. (iii) An interferometer narrowbanded at 
/ = A/27r ~ 1000 Hz with bandwidth A/ = e/27r ~ 50 
Hz (and so denoted "1000"), which has a noise mini- 
mum of ~ 1 X 10^^'*/\/Hz and parameters R — 0.995, 
p = 0.9985, 4) = 1.556. For these three configurations 
the amplification factor is 



.4-wB = 0.27 , A^oo = 1-4 , -4iooo = 5.7 



(4.36) 



We have chosen to compute the signal loss at the opti- 
cal resonance so the signal field u'^ in the arm cavity will 
be on resonance, just as the carrier field is. This allows 
us to translate our carrier-field results over to the signal 
field with only one change: the influence of the tilts of 
the SRM and ITM must be multiplied by the amplifi- 
cation factor A. Therefore, the fraction 5s of the signal 
field's power that is in the tilt-induced parasitic modes 
is 5s — a^^ a\j ~\- CKi5fl). The influence ctj^ of the 

ETM is the same as in the case of the carrier, which we 
have already dealt with, so we shall ignore it here and fo- 
cus on the two mirrors that make up the SR cavity: the 
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ITM and the SRM. If they both have the same tilt an- and 
gles 9 (but about uncorrelated axes) so a\i — a\gj^ = a\, 
then these tilts produce a fractional loss of signal power 
and thence fractional increase in gravitational-wave noise 
given by 



5Sh 
Sh 



(4.37) 



signal 



This is greater by the factor J? than the noise (|4.32() 
due to loss of carrier light into parasitic modes, and cor- 
respondingly to keep this fractional increase of Sh below 
one per cent requires controlling the ITM and SRM tilts 
to an accuracy 1 /A of that in Eq. (|4.33() : 



qMH WB 
^1% 



7 X 10"* 



osph WB 
^1% 



30 X 10"* , 



ef% '™ > 1-4 X 10-* , e^^p^ > 6 X 10-* 

5MH 1000 > 4 ^ ^0-8 



jsphiooo > 10-8 

i /o 



(4.38) 

For the narrowbanded interferometers these limits are 
only lower bounds on Oi% because of the failure of the 
geometric optics limit. As we have seen, the Fresnel 
length for light trapped in the SR cavity is about 1/2 
to 1/3 of the 95-percent-power beam diameter, so trans- 
verse spreading of the light will reduce somewhat the SR 
cavity's amplification factor A and thence the influence 
of tilt on the beam asymmetry. We guess that this re- 
duction might increase 6*^%"^ by a factor of order 2 over 
the geometric-optics limit, Eq. H4.38|l : but since A is only 
about 1 for 6*^% , we guess that there is little increase in 



p500. 

^1% ' 



flMH 500 
qMH 1000 



1.4x10-*, ~ 6 X 10-* , 

0.8 X 10-* , ef^^ ~ 3 X 10-* 



(4.39) 



4. Transmission of Signal Light Through SR Cavity 



df'/d(t) 



4i?p2 



1 



(4.42) 



Rp"^ + 2VRp cos 2(1) 

Here we have used Eq. (11) of Buonanno and Chen 
for f', with the factor e"^^'"^ in the numerator removed, 
so as to obtain the transmissivity that carries the field 
from an input transverse plane to an output transverse 
plane in the presence of the mirror tilt (which gives 
its dependence on transverse position r). For our three 
interferometer configurations, the values of B are 



B- 



WB 



15, B 



500 



33 , Bi 



66 



(4.43) 



When both ITM and SRM are tilted through the same 
angle 9 about uncorrelated axes, the total power trans- 
ferred into parasitic modes is twice as large as Eq. (|4.40|) 
[i.e., 6t is twice (|4.40ll ]. and correspondingly the fractional 
increase of gravitational wave noise (due to both loss of 
signal power and increase of shot noise) is 



Sh 



St=B'k'{r') 



(4.44) 



transmission 



Inserting the above values for B and (r^) and insisting 
that Sh not increase by more than one per cent, we obtain 
the following constraints on the ITM and SRM tilt angles: 



qMH WB 



1.6 X 10-* 
f% > 0.7 X 10-* 

IH 

1% 



nsph WB 



2.4 X 10" 



^MH 1000 > ^ ^0-8 



^sph 500 > ^ ^ ^ 

01% > 0.6 X 10-* 



(4.45) 



For the narrowbanded interferometers, the failure of 
the geometric optics limit dictates that these estimates 
of 6*1% are lower limits; hence the ">" . As in the case of 
signal light in an arm cavity reflecting off the SR cavity, 
so also here, we guess that these estimates are fairly good 
for narrow banding at 500 Hz and are roughly a factor 2 
too severe at 1000 Hz, so 



When the ITM or SRM is tilted through an angle 
0, producing a spatially dependent mirror displacement 
Sz = 9rsinip, it alters the SR cavity's transmissivity by 
a spatially dependent amount 5t' = [df /d(j))kSz, in the 
geometric optics limit. When an undistorted signal beam 
Mo passes through this spatially variable transmissivity, 
a fraction 



{uoASt'/t'I'uo 



1_ 



-B'eir' 



(4.40) 



gets transferred to parasitic modes. Here fc = 27r/Ao is 
the wave number, (r^) = (uo^r'^uo) is the beam's mean 
square radius, which has the values 



.95cm) 



for fiducial mesa beam 



(4.41) 



(4.70cm)'^ for fiducial Gaussian beam. 



/)MH 500 
^1% 

qMH 1000 
^1% 



0.7 X 10-* 
0.7 X 10"* 



qsph 500 
^1% 

qsph 1000 

^1% 



1.1 X 10-* 

i 1.1 X 10-* . 

(4.46) 



Equations (|4.45|) for wideband interferometers and 
KM for narrowband [and for ^"°"] are the 

most severe of all our tilt constraints. 



5. Reference light for baseline readout 

The baseline design for the advanced LIGO interfer- 
ometers includes an output mode cleaner and homodyne 
readout. As we have seen [Eq. H4.31|l ]. the mode cleaner 
makes Sh insensitive to (first-order) losses of reference 
power into parasitic modes. 
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6. Transmission of RF reference light through a power 
recycling cavity without an output mode cleaner 

In LIGO-I interferometers, by contrast with the base- 
hne design of advanced LIGO, there is no output mode 
cleaner, and heterodyne readout is used in place of ho- 
modyne; i.e., the reference light is radio-frequency (RF) 
sidebands, transmitted to the output port via the RF cav- 
ity (which is the same as the power recycling cavity since 
LIGO-I has no SR mirror). In this section we shall com- 
pare our approximate analysis with more careful analyses 
of LIGO-I, by analyzing this type of readout. 

Suppose that the PRM or ITM is tilted through an 
angle 9 and is thereby given the space-dependent dis- 
placement 5z = Orsinip. Then, in the geometric op- 
tics limit, the RF reference light acquires, when passing 
through the PR cavity, a space-dependent phase shift 
(jF/7r)fc2(5z,^^ where T is the cavity finesse. The refer- 
ence light emerging from the cavity therefore has the form 
mJ, — itoe*'-'^/'^^'^''^, for which the fraction of light power 
in parasitic modes is 

uo, (^2^ ' kHz^ - i (^^fc) ' {r^)e^ . (4.47) 

When both the PRM and the ITM are tilted through the 
same angle 9 but around uncorrelated axes, the parasitic 
mode power is twice as large [so Sr is twice (|4.47ll ]. and 
the fractional decrease of signal power, and correspond- 
ing increase of gravity-wave noise, due to the loss of this 
reference-light power, are then 

(f^) =S.= ('-^k)\r^)9^. (4.48) 

\ ^li / reference \ ~ / 

For the LIGO-I interferometers is rather large, T ~ 
105, which produces a strong sensitivity to mirror tilt. 
By contrast, for advanced LIGO J- is smaller, T ~ 50, 
which (as we shall see) compensates for the larger beam. 
The result would be about the same sensitivity to tilt 
as for LIGO-I, if there were no output mode cleaner in 
advanced LIGO and heterodyne readout were used. 

Inserting the LIGO-I finesse ~ 105 and mean square 
beam radius^'^ (r^) = = (2.9cm)^, and constraining 
the increase of Sh to less than one per cent, we obtain 
the following constraint on the ITM and PRM tilts: 

6*1^^ ~ 0.9 X 10"* . (4.49) 

This differs by a factor ~ 1.5 from the result of a much 
more careful computation by Fritschel et. al. 17], which 



gavei4 9];^^°-^ = 1.4 X 10"^ 

If advanced LIGO were to use an RF readout and no 
mode cleaner, the RF cavity would be bounded by the 
PR, SR, and ITM mirrors (with each round trip of light 
encountering the ITM twice). Uncorrelated tilts of all 
three three mirrors would lead to 6Sh/Sh — Sr twice 
as large as expression (|4.48|) . Inserting the advanced 
LIGO finesse ~ 50 and mean-square beam radius 
(r^) = (6.95cm)^ (fiducial mesa) and (4.70cm)^ (fidu- 
cial Gaussian) , and constraining the increase in Sh to no 
more than one per cent, we would then obtain 

9^^ ~ 0.4 X 10-« , 91"^^ ~ 0.6 X 10-« , 

As we have seen, an output mode cleaner will remove 
this increase of Sh, making these constraints no longer 
needed. 



J. Decrease in Signal Strength Due to Mirror 
Figure Errors 

The loss of signal strength and resulting increase in 
Sh, caused by mirror figure errors, is given by the same 
equation Sh 1 + 2(5o + 2(5^ + St + {Sr or 0} as for mirror 
tilt [Eq. H4.31|l ]. but now Se is the fraction of power in 
the parasitic components of mode due to figure errors 
rather than tilt. 



1. Carrier light in arm cavity 

Deformations of the ITM and ETM, with optimized 
tilt compensation, drive a fraction = PicE + l^ici i^to 
parasitic modes [Eqs. H4.22|l and H4.23|l ]. Assuming the 
same peak-to- valley deformations Az in the two mirrors' 
central regions, we obtain for the fractional increase in 
gravitational- wave noise 

(|i)"".2...4,4.0.0l(^)\ (4,50) 

Correspondingly, to keep the noise increase below one per 
cent, we must constrain the ITM and ETM deformations 
to 

AzJ'J/f ~ 15nm (4.51) 

We have not carried out an analysis of the influence of 
the ITM and SRM mirror deformations on the fiducial 
Gaussian arm cavity modes, and so cannot say what the 
analogous constraint is in the Gaussian case. 



In this equation was missing the factor 2, which produced 
fac tor 2 errors in the resulting constraints on tilt; cf. Sec. IV, E, 5 
of [T3 . 

LIGO-I has asymmetric arm cavities; the 2.9 cm that we use is 
the root mean square of the beam radii on the ITM and ETM. 



Their result (end of Sec. 2. A) is 6* < 1.0 X lO"** for the tilts 
in pitch and in yaw, corresponding to a constraint 6 < \/2 X 
IQ—S ^ ]^ 4 X 10~* on the magnitude of the vectorial tilt, for 
an 0.005 fractional decrease in amplitude signal to noise, which 
corresponds to our one per cent increase in 5^. 
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2. Signal light in arm cavity 

As for tilt, so also for figure errors, the SR cavity am- 
plifies the influence of the errors Az by a factor A [Eq. 
(|4.35|) ]. so the ITM and SRM deformations move a frac- 
tion 5s = A^{/3l^j + PIcSr) cavity's signal 
light into parasitic modes. When the two figure errors 
have the same magnitude and we wish to keep the re- 
sulting loss of signal power and increase of Sh below one 
per cent, this gives rise to constraints on the ITM and 
SRM figure errors that are 1 /A more severe than H4.51|l . 
Using the values H4.36|l of A for our three interferometers 
(wide-band, narrowbanded at 500 Hz and narrowbanded 
at 1000 Hz), and increasing the limit for the 1000-Hz 
narrowbanded case by a factor 2 due to failure of the 
geometric-optics limit [cf. Eq. H4.39|l and associated dis- 
cussion], we obtain the constraints 



Azj^f > lOnm , with a guess of ~ lOnm , 
Azj^f 1000 > 2.6nm , with a guess of ~ 5nm . 

(4.52) 



These are the most serious of our constraints on the 
mirror figures of advanced interferometers, and they are 
the same for fiducial mesa and fiducial Gaussian beams, 
because transmission through the SR cavity is governed 
(at least roughly) by geometric optics. The one small 
difference is that the central region over which the peak- 
to-valley deformations are constrained (the region con- 
taining ~ 95 per cent of the light power) is larger for the 
fiducial MH mirrors (central radius about 10 cm) than 
for the fiducial Gaussian mirrors (central radius about 8 
cm) . 

The mirror-figure constraint (|4.54|) is sufficiently se- 
vere, at least in the case of narrowbanded interferome- 
ters, that it might be worth considering reducing the de- 
generacy of the SR cavity by making the entrance faces 
of the ITM's into lenses that bring the signal light (and 
inevitably also the carrier light) to a focus somewhere 
near the SRM (and PRM). Since the constraint 
is the same, whether the mirrors are MH or spherical, 
this recommendation is not dependent on whether mesa 
beams are implemented. 

4- Reference light for baseline readout 



3. Transmission of signal light through SR cavity 



By the same analysis as for mirror tilt (Sec. IIV 1 4p . 
deformations 5z{x,y) of the ITM and SRM by the same 
peak-to-valley amounts Az produce an increase in Sh 
given by 



Sh 



transmission 



St = 2B^k^{{Szf) 

-B^^iAz)^ 
4 



(4.53) 



[cf. Eq. with (((5z)2) = ((6'r sin ^)2) = ^{r^)e^]. 

Here {{Sz)'^) is the mean square deviation of ITM or SRM 
height from the desired figure, and we have approximated 
this by half the squared amplitude of mirror height fluc- 
tuations, which is 1/8 the square of the peak to valley 
height fluctuations, (Az)^/8. Inserting the values of B 
for our three interferometers [Eq. (|4.43l) ] and requiring 
that Sh increase by no more than one per cent, we ob- 
tain the following constraints on the ITM and SRM peak 
to valley deformations: 



A .^MH WB 
A .,MH 500 
A ,MH 1000 



. sph WB 



2nm, 



Az^y^ ^'^^ < Inm, with a guess ~ Inm , 

AzJyJ' ^'^"'^ ^ 0.5nm, with a guess ~ Inm . 

(4.54) 



Here as in Eqs. H4.52|l . (|4.46() and (|4.39() . the breakdown 
of the geometric optics limit in the SR cavity has dic- 
tated a lower limit and a guess for the narrowbanded 
interferometers. 



As for mirror tilt, so for figure errors, the output mode 
cleaner in advanced LIGO makes Sh insensitive to the 
loss of reference light into parasitic modes. 



5. Transmission of RF reference light through power 
recycling cavity without an output mode cleaner 

In the LIGO-I case of no output mode cleaner and 
heterodyne readout with RF sideband light, the loss of 
reference light to parasites does increase Sh- By the same 
argument as for mirror tilt fSec lIVIGl) . the net increase 
of Sh due to deformations 5z{x, y) of the ITM and PRM 
isi5 



5Sh 

/ reference 



2T, 



(Az)' 



(4.55) 



[cf. Eq. J^mi with (((5z)2) = ((6'r sin ^)2) = ^{r^)e^]. 
For advanced LIGO with RF readout and no mode 
cleaner, there is an extra factor 2 in Eq. H4.55(l due to 
change of configuration of the RF cavity (the added SRM 
and the light encountering the ITM twice on each round 
trip). Inserting the finesses of LIGO-I {J- ~ 105) and 
advanced LIGO {J- ~ 50) and insisting that Sh not be 



15 In H the factor 2 in (2J^fc/7r)^ was missing, which produced 
factor 2 e rro rs in the resulting constraints 14.561 on Az; cf. Sec. 
IVF4 of Q. 
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increased by more than one per cent, we obtain the fol- 
lowing constraints on the central-region peak-to-valley 
deformations of the PRM and SRM: 

Az\-^^°-^ = 0.5nm, 

Az^^ = Azl^^^Q.bnm. (4.56) 

These constraints are rather severe, but for advanced 
LIGO the constraint is actually irrelevant because of the 
output mode cleaner (and the planned use of homodyne 
readout rather than heterodyne). 

V. CONCLUSIONS AND FUTURE RESEARCH 

The thermoelastic benefits of MH mirrors are suffi- 
ciently great, and the tightened constraints that they 
place on mirror figures, positions and tilts are sufficiently 
modest, that MH mirrors have been adopted by the LSC 
as an option for advanced LIGO, and will be incorporated 
into future modeling along with spherical mirrors. 

Two methods of manufacturing MH mirrors are be- 
ing explored by the LIGO Laboratory: diamond cutting, 
and evaporative coating. The immediate goal is to de- 
termine the accuracy and reproducibility with which the 
desired MH mirror figures can be produced by each of 
these methods. 

If MH mirrors are to be considered seriously for ad- 
vanced LIGO, it is necessary to develop laboratory ex- 
perience with them and with the mesa-beam optical cav- 
ities that they produce. An experimental effort in this 
direction has been initiated by Phil Willems of the LIGO 
Laboratory at Caltech. 

As we have discussed in Sees. IIV HI and IIV J 31 the 
constraints on mirror figure in the recycling cavities are 
rather worrisome for both spherical mirrors and MH mir- 
rors. More accurate studies of this are needed. As 
one aspect of these studies, it is important to check our 
geometric-optics-based claims that the fiducial MH recy- 
cling cavities are not much more sensitive to mirror figure 
errors than the fiducial Gaussian-beam cavities. 

If the recycling-cavity-induced constraints on mirror 
figure are found to be as serious as our estimates sug- 
gest, then it seems worthwhile to carry out studies of the 
option of converting the input faces of the ITM's into 
lenses that make the recycling cavities much less degener- 
ate. PhD students in Thome's research group may carry 
out these studies as part of the broadening experience 
that is a standard part of their education. 

Bill Kells has suggested the possibility of operating the 
advanced LIGO interferometers initially with spherical 
mirrors and Gaussian beams, and later switching to mesa 
beams by altering only the ETM's. It seems to us that 
this option is worth detailed study. If the ITM input faces 
are turned into lenses that reduce the near degeneracy of 
the recycling cavities, then it might be possible to keep 
the lenses weak enough that the figures of the recycling 
mirrors can be the same for Gaussian and mesa beams. 



while still relaxing the mirror figure constraints to an 
acceptable level. This needs study. 
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APPENDIX: APPROXIMATE FORMULAE FOR 
MESA MODES 

By inserting into Eq. (|2.5|) the asymptotic expansion 
(at large agument) for the Bessel function 

Io{z) = -^e^ , (1) 

setting 6 = 1, To = Z) — ^, and 

a = v/(l + *)/2 , (2) 
we bring Eq. (|2.5(l into the approximate form 

U{D,r) = J"" ^D~^exp[-a^r - D + ■ 

(3) 

By expanding ^/D ~ ^ as a power series in £^/D up to 
some order n, and then performing the integral in Eq. 
(O analytically, we obtain expressions for U{D,r) with 
various accuracies. The least accurate expression, n — 
(obtained by setting vTJ^~J = a/D, integrating, and dis- 
carding a term proportional to erfc(Q!r) that is negligible 
compared to erfc[a(r — D)] at the relevant radii, r ^ D 
or larger) is 

Uo{D,R)^^J-eTic[air-D)]. (4) 
V r 

Here erfc(2:) is the complementary error function, 1 — 
erf(z). At order n = 3 we get a much more accurate 
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FIG. 10: Comparison of the moduli of the exact flat-topped 
mode 17(4, r), the zero-order approximation C/o(4, r) [Eq.^J], 
the third-order approximation f/3(4, r) [Eq. JKJ], and the 
asymptotic approximation to the zero-order approximation, 
Uoaij^yv) [Eq. 0]. Top: linear plot; bottom: logarithmic 
plot. 



expresssion 



X {tt [a^{2r^ - lODr^ + 30L»V + 101?^) + a{3r 
X erfc[a(r — D)] — erfc[ar] 

+2^/¥ 



5^)] 



[a^ (11 ~ 4:Dr + r^) + 1] } 



(5) 



Here erfc is the complementary error function, erfc(z) = 
1 - erf(z). 

The two approximations C/o(4, r) and 1/3(4, r) are com- 
pared with the exact mode C/(4,r) in Fig. ^] At r > 2, 
Us is highly accurate; at r > 5, Uq is highly accurate. The 
analytic formula Q for Uq shows that the diffraction- 
induced tail of this mode falls off very rapidly outside 
r = D — slightly more rapidly than would the tail of a 
Gaussian centered on r = D. Of course, this is to be ex- 
pected since U{D,r) is constructed from a superposition 
of Gaussians that are centered on radii < D. 

A useful but cruder approximation to U{D, r) can 
be obtained by inserting the asymptotic expansion of 
erfc[Q;(r — D)] into Eq. |@J for Uq. The result is 



UoaiD,r) 



which has a modulus 



(6) 



\Uoa{D,r)\ 



V2r 



r-D 



(7) 



This modulus is accurate to within several tens of per 
cent in the regime of interest for diffraction, D ^ A and 
i? > 6 (in units of h). By inserting this approximation 
into Eqs. (|2.11|) and H2.7|l for diffraction losses in the 
clipping approximation, we deduce that the diffraction 
losses for mesa beams scale with beam-spot radius D as 



-2RD/b'^ 



(8) 



Here R is the mirror radius and we have restored the 
diffraction lengthscale b, which was set to one throughout 
this appendix. 
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